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We use a number of large- N limits to explore the competition between ground states of square lat- 
tice doped antiferromagnets which break electromagnetic U(l), time-reversal, or square lattice space 
group symmetries. Among the states we find are d-, (s* + id)-, and (d x 2_ y 2 + id xy )-w&ve supercon- 
ductors, Wigner crystals, Wigner crystals of hole pairs, orbital antiferromagnets (or staggered-flux 
states), and states with spin-Peierls and bond-centered charge stripe order. In the vicinity of second- 
order quantum phase transitions between the states, we go beyond the large- N limit by identifying 
the universal quantum field theories for the critical points, and computing the finite temperature, 
quantum-critical damping of fermion spectral functions. We identify candidate critical points for the 
recently observed quantum-critical behavior in photoemission experiments on Bi2Sr2CaCu2 0s+a by 
Valla et al. (Science 285, 2110 (1999)). These involve onset of a charge density wave, or of broken 
time-reversal symmetry with d x 2_ y 2 + id xy or s* + id pairing, in a d-wave superconductor. It is 
not required (although it is allowed) that the stable state in the doped cuprates to be anything 
other than the d-wave superconductor — the other states need only be stable nearby in parameter 
space. At finite temperatures, fluctuations associated with these nearby states lead to the observed 
fermion damping in the vicinity of the nodal points in the Brillouin zone. The cases with broken 
time-reversal symmetry are appealing because the order parameter is not required to satisfy any 
special commensurability conditions. The observed absence of inelastic damping of quasiparticles 
with momenta (tt, k), (k, n) (with < k < tt) also appears very naturally for the case of fluctuations 
to d x 2_ y 2 + id xy order. 
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superconductors is to assume that all relevant ground 
states can be completely characterized by the manner in 
which they break the global symmetries of the underlying 
Hamiltonian. Related ideas have been discussed by oth- 
ers in Refs. |^-^; for a review by one of us see Ref. [|, and 
for early work we shall extend in this paper see Ref. §. 
The global symmetries are: 

(i) S - the electromagnetic U(l) symmetry, which is bro- 
ken by the appearance of superconducting order; 

(ii) M - the SU(2) spin rotation invariance symmetry, 
which is broken by magnetically ordered states like the 
Neel state; 

(Hi) C - the space group of the square lattice, which we 
will consider broken when an observable invariant under 
S and M, like site or bond-charge density, is not invari- 
ant under space group transformations; and 
(iv) T - time-reversal symmetry. 

Even in this limited approach, a surprisingly rich num- 
ber of phases and phase transitions are possible. While 
at low enough temperatures, every phase is amenable to 
a conventional quasiparticle-like description, anomalous 
behavior can appear in the vicinity of second-order quan- 
tum critical points between the phases. 

This paper will study the competition between phases 
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which break one or more of iS, C, and T symmetries, and 
describe the universal theories of the associated quan- 
tum critical points. We will do this in the context of 
a number of large- N computations which, by construc- 
tion, only produce ground states in which A4 symmetry 
is preserved. Transitions at which A4 symmetry is bro- 
ken are quite important for certain aspects of the physics 
of the cuprates (as we shall discuss later in this section), 
and the absence of explicit computations for where such 
transitions may occur is the primary limitation of our 
approach. 

A brief discussion of mainly the large-iV results has 
appeared in an earlier papem. 

A global perspective on our approach is provided by 
the schematic phase diagram in Fig. |l| (see also the dis- 
cussion by ZaanenEl on related phases) . This sketches the 
qualitative evolution of the physics as a function of the 
hole doping, <5, and N, the size of the symmetry group of 
spin rotations, M: we will primarily consider models in 
which this symmetry group is generalized from SU(2) to 
Sp(2A^) (specific details in Section |J). 

Let us first discuss the physics at S = as a function 
of N. For small N, including the physical case, N = 1, 
we know that the ground state has magnetic Neel order, 
and so breaks M. symmetry. This symmetry is restored 
by a continuous quantum phase transition at the point X. 
Above X, it is believed that one enters a ttataipagnetic 
phase which generically has 'Peierls' ordenTBEO'EJ: in 
such a state all sites are equivalent, but the charge and 
energy densities on the bondsEj have the modulation indi- 
cated schematically by the pattern in Fig. |l|. It is evident 
that such a state breaks only a C symmetry. So the M. 
and C symmetries vanish at the common point X, a phe- 
nomenon not generically expected in a Landau theory- 
like approach suitable for high dimensions, but possible 
in the present low=dimensional system with strong quan- 
tum fluctuations^. Further above the i-pnkit X, other 
phases like the 'orbital antiferromagnetoTS (to be de- 
scribed shortly) are also possible in certain models and 
we will also discuss these. 

The primary purpose of this paper is to study the 
physics for S > 0. Ideally, we should do this along the 
path A2 in Fig. [|, which meets the 5 = line below 
X. Instead, we will offer a controlled, quantitative the- 
ory along the path A\ which meets the 8 = line above 
X. The implicit assumption underlying such a strategy is 
that the quantum critical point X is "close" to the regime 
of physically relevant parameters. Then we can expect 
that the phases with M. symmetry preserved along the 
path A\ are related to the phases with M. symmetry pre- 
served accessed by increasing 5 along the path A2 . 

This paper describes the intricate interplay between 
C, S, and T symmetry breaking along A%: the body of 
the paper contains a large number of phase diagrams as a 
function of S and a dimensionless measure of the strength 
of the long-range Coulomb interactions, and all of these 
lie along A\. Over a significant regime of parameters, we 
find that C symmetry is broken at smaller values of 6; 
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FIG. 1. Schematic phase diagram of doped antiferromag- 
nets. See Section |j| for further general discussion and Sec- 
tion |n] for specific quantitative phase diagrams along the line 
A\. The vertical axis represents the size, N, of the symme- 
try group of spin rotations, M. Although this parameter is 
not experimentally variable, we propose that a similar phase 
diagram would be obtained as a function of the ratio of the 
second (^'j]^|%^1Hi>-0 neighbor exchange interactions. There 
is evidenced TJCnij'E-j that the Peierls order shown above X at 
8 = is also found in J' — J antiferromagnets. The magnetic 
M symmetry is broken in the hatched region, while C sym- 
metry is broken in the shaded region; there are numerous ad- 
ditional phase transitions at which the detailed nature of the 
A4 or C symmetry breaking changes - these are not shown. 
For 8 = 0, M symmetry is broken only below the critical 
point X, while C symmetry is broken only above X. Over a 
significant parameter regime, and for not too small 8, the C 
symmetry breaking appears in the stripe patterns shown, with 
accompanying anisotropic superconductivity which breaks 5. 
For some other parameter regimes (as in Fig [ll]) , the C sym- 
metry breaking is realized by orbital antiferromagnetic (or 
staggered flux) order: such C symmetry breaking cannot sur- 
vive all the way down to the point X. The superconductivity 
is pure d-wave only in the large 8 region where C and M 
are not broken. The smaller 8 region of the superconductor 
which preserves C and M can also exhibit (d x 2_ y 2 +id xy )- or 
(s* + id)-wave superconductivity. 

for very small values of 6, the C-broken phase is an insu- 
lating Wigner crystal-like staia, but for larger S we-pb-. 
tain a state with co-existingaO stripe charge orderE3~E3 
and superconductivity, as sketched in Fig. |l|. Moving 
to smaller values of N in Fig. |l|, we expect a transi- 
tion to a state with A4 symmetry broken which is not 
contained in our computations here. The magnetic or- 
der appears in a background of charge stripe order that 
is present on both sides of the transition, and the spin 
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polarization is therefore expected to be incommensurate 
and collinearn. It is important to contrast this magnetic 
order from the incommensairate "spiral" states that were 
considered some time agcc3: such states have coplanar 
spins and do not require co-existing charge order. In 
contrast, the incommensurate, collinear spin states must 
coexist with charge stripe ordercil, and there is evidence 
that the spin incommensttration observed experimentally 
is indeed of this typeE3~E3. Stripe charge order has been 
discussed first in^he-context of mean-field theories for the 
Hubbard modeEHcS, a brief comparison of the se resu lts 
with the present theory will be given in in Sec. |EI A3| . 

The order parameter for the collinear s pia or dering, 
discussed above, is an ordinary 0(3) vectorESE§E3, iden- 
tical to that for the transition at the point X at 5 = 0. To 
understand the transition at which M. symmetry is re- 
stored at S > 0, we need to explore the possibility of the 
magnetic order parameter coupling to gapless fermion 
excitations. In all our charge stripe states, we find a 
strong pairing tendency between the holes, and as a con- 
sequence, the fermion excitation spectrum is either fully 
gapped, or has gapless excitations only at special points 
in the Brillouin zone. Assuming that momentum con- 
servation prohibits the coupling between the magnetic 
order parameter and the gapless fermion excitations (if 
present) , we arrive at the conclusion that the spin disor- 
dering transition at S > is inprecisely the same univer- 
sality class as that at X. Fig. |l| contains a line, emerging 
from the point X, along which a transition to M. sym- 
metry restoration takes place. The gist of our arguments 
above is that the universality class of the transition all 
along this line is likely to be identical to that at X . The 
implications—for such a scenario for experiments (espe- 
cially NMREll) has been reviewed recently in Ref. [|. A 
related scenario, and quantitative comparisons withrex- 
periments, has been provided recently by Morr et aZ.Eil. 

Returning to the physics along At, we briefly catalog 
the properties of the states found. Further details ap- 
pear in Section 0, but the reader is urged to glance at 
the phase diagrams in Figs |^-|TT| at this stage. It is also 
worth noting explicitly here that all of these phase dia- 
grams were obtained in the large N limit, and the precise 
numerical values of the parameters at the phase bound- 
aries are not expected to be accurate for SU(2): never- 
theless, the general topology of the phase diagrams, and 
the trends in stability between the various phases, are 
expected to be realistic. 

(i) Superconductors with C symmetry: These appear for 
large 5, and for a large region of parameters the Cooper 
pairs-are in a <i-wave state. We also fi nd a (s* + ic?)-wave 
state: 



a (d x 2_ y 2 + id xy )-w&ve stat 



both of which 

break T symmetry, but preserve C symmetry; the. latter 
state has a non-vanishing spin Hall conduct anceE3. 
(ii) Spin-Peierls: An insulating spin-Peierls state at S — 
was discussed above: it breaks only C symmetry. We 
also obtainedtj for a range of 5 > a superconducting 
state with precisely the same pattern of C symmetry 
breaking; naturally, S is also broken in such a state. The 



signature of both states in neutron scattering would be 
the same: all sites are equivalent, but there is a modu- 
lation in the energy and charge densities on the bonds 
with a period of 2 lattice spacings. 

(Hi) Stripes: The striped states are similar to the su- 
perconducting spin-Peierls states above, but all sites are 
no longer equivalent. The states have apx 1 unit cell 
and the holes are concentrated on a strip of width q; 
both p and q were always found to be even. The width 
q regions form strong one-dimensional superconductors 
(Luther-Emery liquids), and coupling between these hole- 
rich regions leads to anisotropic superconductivity. The 
hole pairing also, always prefers stripes in a bond-centered 
configuratior£3£l: the ground state possesses a reflection 
symmetry about the centers of certain columns on-bonds, 
but not about any^himn of sites. Distinguishingcl bond- 
and site-centeringE3~E£l is an important issue to be re- 
solved by future experiments. 

(iv) Wigner crystal: An insulating state for 5 > is the 
familiar Wigner crystal of holes, which breaks only C 
symmetry. This state appears when the strength of the 
Coulomb interactions is large enough. 

(v) Pair crystal: An alternative insulating state can ap- 
pear at small S > and weaker Coulomb interactions. 
The exchange interactions induce pairing of holes, and 
the resulting composites then form a Wigner crystal to 
minimize the Coulomb repulsion. The underlying square 
lattice can induce strong distortions on the usual triangu- 
lar structure of the Wigner crystal, so that the state can 
look like a striped configuration with an additional lon- 
gitudinal charge modulation along each hole-rich stripe. 

(vi) Orbital antiferromaanet: This also known as the 
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"staggered flux" statetTlia. There are staggered, circu- 
lating currents around each square lattice plaquette, and 
the state breaks T symmetry. The unit cell has two sites 
and so C symmetry is also broken, but a combination of 
translation by one site and time-reversal remains unbro- 
ken. At half-filling, this state has gapless fermionic ex- 
citations at nodal points along the diagonals of the Bril- 
louin zone, like the d-wave superconductor. The state 
with co-existing orbital antiferromagnetism and d-wave 
superconductivity has broken C symmetry (the unit cell 
has 2 sites) and is distinct from the (s* + id)-wave or 
(d x 2_ y 2 + id xy )-v?&ve superconductor. 

The competition between C, M. and S breaking lead- 
ing to phases with co-existing orders-/ as in (ii) and (Hi) 
above) has been discussed by Zaanenld on phenomcnolog- 
ical grounds. However, he focuses mainly on the bosonic 
order parameters, while fermionic excitations will play 
an important role in our considerations. 

Phases related to those in (i), (Hi), (iv) and (v), and 
associated phase diagrams^, appear in the work of Kivel- 
son, Fradkin and Emerycj. These authors use liquid- 
crystal like pictures, in which quasi-independent one- 
dimensional Luttinger liquids are allowed to fluctuate 
transversely, and obtain qualitative phase diagrams. In 
contrast, our work is intrinsically two-dimensional, and 
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the quantum fluctuations are tied more strongly to the 
underlying lattice sites; further we shall obtain quantita- 
tive results for phase diagrams, albeit in a large N limit. 

It is interesting that none of the states above is an or- 
dinary Fermi liquid: such a state appears to be always 
unstable to some ordering induced by breaking one or 
more of C, S, M., or T symmetries. In this respect o 
results are similar to recent results of Honerkamp et al. 
and Ledermann et alSB; however, they find renormal- 
ization group instabilities of the Fermi liquid to states 
somewhat different from those discussed above. 

We conclude our discussion of mean-field theory by 
noting that a separate study of C symmetry breaking in 
doped antiferropagnets has been carried out recently by 
Stojkovic et alEB: they examined the competition be- 
tween stripe and Wigner crystal-like phases in a semi- 
classical theory of hole dynamics. 



A. Quantum phase transitions 

We have already discussed above the nature of the M. 
ordering transition in-Kig. |l[ and its possible relationship 
to NMR experimentsoEl Here we shall explore the na- 
ture of the quantum phase transitions between the phases 
found along A\ (which do not involve order parameters 
associated with A4 symmetry), and their possible rela- 
tionship to quantum critical scaling observed in a recent 
photoemission experimented. 

Many of the transitions in our phase diagrams are first 
order. These do not have interesting fluctuation spectra, 
and we will not consider them further. We will consider 
second order transitions at <5 > 0, in which the ground 
state is superconducting on both sides of the transitions 
(see Fig ||) . One of these states is a d-wave superconduc- 
tor, while the other is denoted as superconducting state 
X in Fig|^; we will discuss different candidates for state X 
below. Such second order transitions fall into two classes 
depending on the behavior of the fermion spectra in the 
vicinity of the nodal points of the cZ-wave superconduc- 
tor. These points are afl (±K, ±K) with K = 0.391tt 
(at optimal doping), and throughout the remainder of 
this paper, unless noted otherwise, we will be implicitly 
referring to the fermionic Bogoliubov quasiparticles in 
the vicinity of these points. (We will discuss the prop- 
erties of the gapped fermionic quasiparticles near (0, 7r), 
(it, 0) at the end of Section |.) The two classes are: 

(A) There is efficient scattering and damping of the nodal 
fermionic quasiparticles, and as a result the fermionic 
spectral function obeys 'naive' quantum critical scaling 
(sec ([hi]) below), of the type observed experimentallyc3. 

(B) the gapless, Bogoliubov, fermionic quasiparticles can 
be neglected in the scaling limit of the critical theory, and 
so their damping appears only upon considering correc- 
tions to scaling, and vanishes with a super-linear powers 
of temperature (T) as T — > 0. 

The simplest of the transitions in class A are those that 
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FIG. 2. Finite temperature (T) phase diagram in the 
vicinity of a second order quantum phase transition from a 
d-wave superconductor as a function of some parameter in 
the Hamiltonian, so (which is possibly, but not necessarily, the 
hole concentration 8). Superconductivity is present at tem- 
peratures below T c , and the superfluid density is non-zero on 
both sides of sq c . The state X is characterized by some other 
order parameter (in addition to superconductivity) which van- 
ishes above a temperature Tx- We will consider a number of 
possibilities for the state X in this paper, including broken 
T symmetry in a (s* + id)-wave or a (d x 2_ y 2 + id xy )-wme 
superconductor, or in an orbital antiferromagnet, and states 
with broken C symmetry with charge density wave order. Of 
particular interest will be transitions for which the scaling 
form ( [l.l| ) applies to the nodal quasiparticles in the quan- 



tum critical region T > T qc . As we discuss in Section [A 



this scaling could continue to apply even above T c provided 
the thermal length associated with proximity to the quantum 
phase transition at s — so c remains smaller than the phase 
coherence length. It is the proposal of this paper that the 
high-temperature superconductor studied in the experiments 
of Ref W(j is in the vicinity of so c - Indeed, this system could 
have so > so c , so that the ultimate ground state is an ordinary 
d-wave superconductor — the fermion spectrum then exhibits 
consequences of fluctuations into state X at T > T qc . 



involve time-reversal symmetry breaking in the <i-wave 
superconductor: transitions from a <i-wave superconduc- 
tor to a state X which is either a (s* + icQ-wavecaEl or 
a (d x 2_ y 2 + id xy )-wa.v£jB3 superconductor. It is impor- 
tant to note that both these transitions occur for only 
at a finite attractive coupling in the s* or d xy pairing 
channels. This is to be contrasted to pairing instabili- 
ties of a Fermi liquid, which would occur at infinitesimal 
attraction in either channel. However, when the parent 
state is a 4^2 _ y 2 -wave superconductor, the vanishing den- 
sity of states at the Fermi level removes the usual BCS 
log divergence in the Cooper pair propagator, and a fi- 
nite attraction is required for further pairing in the s* or 
d xy channels. This finite coupling instability is directly 
responsible for a non-trivial quantum critical point, with 
strong thermal and quantum fluctuations leading to class 
A behavior, whose effects we shall describe and exploit 
in this paper. Much of the more recent discussion (see 
Ref. ^l] and references therein) of (d x 2_ y 2 + id xy )-wave 
superconductivity has focused on the case where this or- 
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der is induced by external perturbations like an applied 
magnetic field, and is motivated partly by the experi- 
ments of Ref. In contrast, in our paper, we are inter- 
ested in the spontaneous appearance of such order, and 
this is necessarily associated with a sharp transition and 
T symmetry breaking. There is a divergent susceptibil- 
ity associated with this symmetry breaking, and this will 
lead to a low energy amplitude fluctuation mode distinct 
from that discussed in Ref. [H]. 

Another transition involving breaking of time-reversal 
symmetry is that between a d-wave superconductor and 
a state X in which d-wave superconductivity and orbital 
antifcrromagnetism coexist. Unlike the above, this tran- 
sition will be shown to be in class B. 

A slightly more complicated transition in class A is 
one in which X involves the onset of C symmetry break- 
ing in a d-wave superconductor: such a transition is the 
boundary of the C-broken region in Fig. []]. However, a 
special condition is required for such a transition to be in 
class A: the charge-ordering wavevector should precisely 
equal the wavevector between two nodal points in the d- 
wave superconductor; otherwise the transition is in class 
B. The theory for such a transition is closely related to 
models for the onset of antiferromagnetism in a d-wave 
superconductor considered recently by Balents et als3. 
We note that others have also discussed quantum phase 
transitions involving stripe or charge density W^ijie order 
in the cuprate superconductors in recent yearad'O; how- 
ever, in contrast to us, these works have either ignored 
interplay with the dynamic properties of the fermionai , 
or focused on transitions in a Fermi liquid ground stateE3, 
not a d-wave superconductor. 

The hallmark of the class A_transition is that in its T > 
"quantum-critical" regior£3 (see Fig ||), the fermion 
Green's function near one of the nodal points of the d- 
wave superconductor obeys 



G f (k,w) 



vk 



k B T {k B T)V* 



(1.1) 



here z is the dynamic critical exponent (z = 1 for the 
specific models solved in this paper) , k measures the dis- 
tance from one of the nodal points of the d-wave super- 
conductor, lu is a measuring frequency, v is a velocity 
(for z = 1), A is an overall amplitude, rjf is a universal 
anomalous dimension, and is a universal scaling func- 
tion of its two arguments. We emphasize that unless the 
system happens to be precisely at the quantum-critical 
point, which is genetically not expected to be the case, 
the scaling form (LI) will eventually fail as T — > 0. In- 
deed, if define a lower crossover temperature T qc so that 
Q holds for T > T qc , then T qc ~ \s ~ s 0c \ zv where s 
is some coupling constant in the Hamiltonian, the T = 
quantum critical point is at sq = sq c , and v is the usual 
correlation length critical exponent. For T < T qc , nor- 
mal Bogoliubov quasinacticle behavior emerges, and this 
is indeed observedcJcJCJ to be the case experimentally 
at very low T. We emphasize further that it is not even 



necessary that the point sq = sq c be in an experimentally 
accessible parameter regime. So if we are considering the 
d-wave to (d x 2_ y 2 + id xy )-~w&ve transition (for definite- 
ness), then it is not required (although it is permissible) 
that the true ground state of a cuprate compound be a 
(d x 2_ y 2 + id xy )-wetve superconductor over some doping 
regime. It is only necessary that a (d x 2_ y 2 + id xy )-wetve 
superconductor be close enough to the physical regime, so 
that dynamic fluctuations to (d x 2_ y 2 + id xy )-w&ve order 
are apparent in the quantum-critical regime. Postulating 
the existence of a (d x 2_ y 2 +id xy )-w&ve ground state some- 
where in parameter space is then a powerful theoretical 
tool for obtaining a controlled description of this interme- 
diate temperature regime. When this paper was almost 
complete, new experimental evidence for broken T sym- 
metry near defects was reportedEJ; these results support 
the hypothesis that the bulk energy of a superconductor 
with broken T symmetry is not ..very much higher than 
that of a d-wave superconductor^, and that a quantum 
phase transition between these states may indeed be near 
the experimentally accessible parameter space. 

One of the purposes of this paper is to develop a 



method for computing the universal function in (1.1) 
for the various transitions in class A noted above. We will 
find that two different methods are necessary, depending 
on the frequency/ wavevector regime being accessed. For 
vk 3> k B T or Tilo 3> k B T, a straightforward resummation 
of a renormalized perturbative expansion suffices. How- 
ever, for vk < k B T and Tiuj < k B T, an entirely different 
approach has to be developed. Now there is strong damp- 
ing induced by scattering between thermally induced ex- 
citations, and we compute it in a self-consistent theory 
of excitations scattering via a renormalized, temperature 
dependent T-matrix. 

It is important to keep in mind that the class A tran- 
sitions being considered here have long-range supercon- 
ducting order on both sides of the quantum critical point. 
The order parameter associated with the transition in- 
volves either C or T symmetry breaking, and has no di- 
rect relationship to S symmetry. The main role of the 
superconducting order is to define the bare spectrum of 
the fermion excitations which then interact with the crit- 
ical order parameter fluctuations. In the experimentsCj, 
scaling related to (1.1) is also observed above the su- 



perconducting transition temperature, T c (which is quite 
distinct from T qc and could be either above or below it); 
indeed there is no signature of T c in the photoemission 
spectrum (while below T qc there is a crossover to con- 
ventional quasiparticle behavior). Our quantum-critical 
theory entirely neglects the fluctuations of the supercon- 
ducting order itself, but this does not limit its applicabil- 
ity to below T c ; rather, we only need to impose the more 
limited constraint that the phase coherence length of the 
superconducting order parameter is larger than the in- 
elastic scattering length of the class A transition with C 
or T symmetry breaking. This constraint is automati- 
cally satisfied below T c , and can easily be satisfied over 
a wide range of temperatures above T c . Indeed, the lat- 



5 



ter length, by (1.1), decreases as ~ l/T 1 / 2 , and so the 
constraint becomes less stringent as T increases. 

Finally, we discuss the important issue of the gapped 
quasiparticles well away from the nodal points of the d- 
wave superconductor. The fermionic quasiparticles with 
momenta (it, k), (k,Tr) (with < k < tt) have a p<?p-| 
zero excitation energy which is a minimum pUcarcS'Ca 
k w 0.187T at optimal doping. ExperimentsCjO clearly 
indicate that these quasiparticles are sharply defined at 
all temperatures below T c , and do not show any sign 
of "quantum-critical" damping (we thank M. Norman 
and M. Randeria for emphasizing this to us). So the 
order parameter fluctuations discussed above, which are 
responsible for the quantum-critical damping near the 
nodal points (±2f , ±K) , clearly cannot couple efficiently 
to the quasiparticles on the lines between (tt, tt) and 
(71-, 0), (0, tt). For the case of a transition from d-wave 
to (d x 2 ,,2 + id xv )-wsive superconductivity (discussed in 



Section III A ) this is, in fact, very naturally the case: the 
d xy order parameter ~ sin k x sin k y vanishes when either 
:. For the transition from d-wave to 



tt or k 



(s* + id)-wave (also discussed in Section III A), the s* 
order parameter ~ (cosk x + cosk y ), and this vanishes 
at the points (tt,0), (0, 7r); however there will be some 
residual coupling as one moves away from these points to 
(tt, k), (k,ir) with k w 0.187T. Finally, for the transitions 
in class A involving the onset of C symmetry breaking (to 
be discussed in Section MB ), momentum conservation 
makes the coupling of the order parameter to fermions 
near (tt,0), (0, 7t) very ineffective: the order parameter 
scatters the fermions to a region of the Brillouin zone 
where the quasiparticles have an even higher energy. 

The outline of the remainder of this paper is as fol- 
lows. In Section || we present the results of the large- Af 
study along A\. The universal theories of the second- 
order quantum phase transitions appear in Section III. 
A summary of our results and a discussion of experimen- 
tal issues is in Section IV. A calculation of the fermion 



damping in a naive renormalized perturbation theory, 
and its failure in the low frequency regime Tiuj < ksT 
is discussed in Appendix [a|. Readers not interested in 
specific details of our results can glance at Figs p[]l~T| and 
move ahead to Section fv[ 



II. Sp(2A0 t-J MODEL IN THE LARGE- iV LIMIT 

For a microscopic investigation of the ground states 
of doped antiferromagnets we start from the usual t — J 
model Ht,j on the sites i of a square lattice, which is 
complemented by a Coulomb Tiv interaction between the 
electrons, Ti = Ti t j + Tiv, 



H t j = [-Uj4 a c 3a + H.c. + J K 

i>j 

Ti V = ^ V ij n i n j- 



(Si 



The electron operators exclude double occupancies 
which expresses the (infinite) electronic on-site repulsion. 
rii = c\ a Ci a is the charge density at site i. We will pri- 
marily concerned with the case where the fermion hop- 
ping, tij , and exchange, Jy , act only when i, j are nearest 
neighbors, in which case fcy = t and Jy = J; however, 
we will occasionally refer to cases with second neighbor 
hopping (t 1 ) or exchange (J'). For the off-site Coulomb 
repulsion Vy we assume weak or no screening since the 
zero-temperature ground states will be either insulating 
or superconducting; in both cases the density of states 
at the Fermi level vanishes. Therefore we will use a 1/R 
decay of the interaction, Vij = V/|Rj — Rj|, and the 
strength of the repulsion is parameterized by V. The 



Vij arc included to count 
J model! 



phase separation ten- 
and play a key role in 



(2.1) 



dency of the t 
our analysis. 

We shall be interested in describing the ground state of 
Ti as a function of its couplings and the average d(j 
concentration S. We generalize the spin symmetry 
from SU(2) to Sp(2A^) and examine the limit of large de- 
generacy N. In the large- N approach the ground state of 
the system can be found in a saddle-point approximation 
(which becomes exact for N = oo). Large- N expansions 
have been applied to a large number of antiferromagnetic 
spin systems as well as to models with doping. The mo- 
tivation for using the symplectic Sp(2A^) generalization 
instead of the more common SU(A r ) variant is that it 
does not rely on a two-sublattice structure of the under- 
lying antiferromagnet and is therefore more appropriate 
for systems where frustration may play a roleo. Further- 
more, the Sp(2A^) approach includes naturally pairing 
of spins which leads to superconducting ground states 
when the system is doped. (Note that both Sp(2A^) and 
SU(2Af) reduce to the usual SU(2) symmetry group for 
N = 1.) Apart from studies of t-J models like ( Oj), 
Hubbaijd-models with Sp(2A^) symmetry have also been 
studiedBEa. 

The behavior of the system in the large- limit de- 
pends on the representation of Sp(2A^) used for the 
spin operators. In the context of SU(A f ) approaches to 
the 2D Heisenberg model especially totally symmetric 
(bosonic) and totally antisymmetric (fermionic) repre- 
sentations have been applied, see e.g. Refs. [l^Jlq^ , |60[ 
For the present problem of the t — J model it turns out 
that a simple large-A^ limit (leading to a saddle point 
in the free-energy functional) exists only for a fermionic 
representation of the spin degrees of freedom. As dis- 
cussed in the introduction, in this case the ground state 
of the undoped (Heisenberg) model does not break spin 
rotation symmetry M.. Instead, translation and rotation 
symmetries C are broken and the state has been shown 
to be a paramagnetic spin-Peierls state which can also 
be considered., as a bond- centered charge density wave. 
Recent workE3 has shown strong evidence for this order 
in the frustrated SU(2) quantum antiferromagnet on a 
square lattice (with J' > 0). 
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Let us now describe the details of the large-iV ap- 
proach. In the following we consider spins transforming 
under the antisymmetric product of m fundamentals of 
Sp(27V), the large-iV limit is taken with m/N constant. 
The spins are represented by fermions f a , a — 1 . . . 2N, 
which transform under the fundamental of Sp(2N). The 
holes are described by spinless bosons b, cf = ffb\. The 
local constraint of the t — J model acquires the form 



flft + bth = m. 



(2.2) 



Here we will only discuss states being half-filled at zero 
doping, m — N . The average hole concentration S is 
determined by 



1 



£/JU° = N(1-S) 



(2.3) 



where N s denotes the (infinite) number of lattice sites. 
Within the Sp(2iV) generalization of the system, the spin 
operators Si become fermion bilinears tim es th e traceless 
generators of Sp(27V), the Hamiltonian (2.1) takes the 
form 



** = £ 
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where we have split Htj = Tit + Hj for convenience. 
The Lagrange multipliers A^ enforce the local occupation 
constraint, and // fixes the average hole density. J a ^ 
denotes the antisymmetric Sp(2iV) tensor: 
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We re mind the reader that Ti given in eq. fl2.4| ) reduces 
to eq. ( |2.l| ) for the "physical" case of N = 1. 

In the limit N = oo at zero temperature the bosons 
bi condense, (fej) = y/Nbi, so bf is the hole density and 
N{\ — bf) = (rii) the charge density at site i. The exclu- 
sion of double occupancies in the N — 1 case is now rep- 
resented by the on-site constraint (m) < N since b\ > 0. 
The long-range nature of Coulomb interaction requires 
the introduction of a background charge of magnitude 



5 on each lattice site for the total system to be charge- 
neutral. The interaction is decoupled by the introduction 
of the link fields Qij with the saddle-point values 



1 



NQa = (J a0 ftfl) = ^{J a,3 4 a c] p ) ■ (2.6) 

1 3 

The Qij represent the complex bond pairing amplitudes. 
We note that Qij = Qji, and the phases of the Qij are 
only fixed up to a global gauge transformation, fa — > 



fie , which leads to Qi 



-2i0 



However, the 



plaquette operator Hq, which can be defined as Hq — 
Q12Q23Q34Q41 f° r tne s ^es 1...4 at the corners of a unit 
square, is a gauge-invariant object. 

At the saddle point the Hamiltonian takes the form 



H t = 



£[-*eW<U° + H.c. 



£ Ai (flfi + Nb " - N ) + M^£ (b 2 i - S) , 



Hj = Y,\- J -f{j aP flfUQ%+^ 
%>j L 

n v ^j2 Nv ^( b "- s ^- 6 y 



N\Q t 



(2.7) 



which is bilinear in the fermions and can be solved by a 
Bogoliubov transformation. The saddle point solution is 
found by minimizing the total free energy with respect to 
bi and Qij at fixed average fermion occupation and hole 
density. The saddle-point equations have been solved nu- 
merically with unit cell sizes up to 32 sites. The Coulomb 
repulsion term in the large- N limit is purely "classical" , 
i.e., it does not involve the fermions. Using the lattice 
Fourier transforms Vk and rik of interaction and charge 
distribution respectively, the Coulomb contribution to 
the energy can be re-written as 



N s NY,V k n k n- 



(2.8) 



The interaction Vk = V J2r^o elkR /|R-l behaves as ~ 
k^ 1 for small k = |k|, but becomes negative around the 
center of the Brillouin zone. (Note that there is no on- 
site contribution from V, Vu = 0, since the occupation 
constraint is already taken care of by the local chemical 
potential A,.) 

The numerical determination of the minimum energy 
configuration in the large- N limit consists of two nested 
loops: (i) Starting from an initial guess for Qij, bi, Xi 
the fermionic Hamiltonian is diagonalized using a dis- 
crete momentum grid (32 2 is sufficient for unit cells 
up to 8 sites). The expectation values (fj a fjp) pro- 
vide new values for the Qij, new bi are obtained from 
&| = 1 — (// Q /f) /N , and the average chemical is ad- 
justed (by a simple bisection step) to match the doping 
level. This is repeated until convergence is reached, (ii) 
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The optimum charge distribution within the unit cell of 
N c sites is found by minimizing the total energy [ob- 
tained in loop (i)] w.r.t. the differences in the chemical 
potentials A,; = A; — Ai. This is a minimum search in 
a (N c — l)-dimensional space and can be performed by 
standard methods. To account for the possible existence 
of more than one saddle point the initial link field values 
Qij are chosen randomly and several sets of initial condi- 
tions are used to identify the saddle point corresponding 
to the global minimum energy. 



A. Ground states at N = oo 

The results of our large- N calculation can be summa- 
rized as follows: First, at 6 — along A\ we find the fully 
dimerized, insulating spin-Peierls (or 2x1 bond charge 
density wave) solution™ in which \Qij\ is non-zero only 
on the bonds shown in Fig. [I]. 

At non-zero doping 6 the "bare" large- N t — J model 
shows phase separation for a large range of parameters 
t/J, see Ref. |^. This tendency to phase separation (into 
a hole-rich region and a fully dimerized half-filled region) 
at V — is an important ingredient of our study. With 
the inclusion of Hy the phase separation becomes "frus- 
trated", as emphasized in Refs. |52"|-f54], and the compe- 
tition of the energy scales t, J, and V leads to various 
kinds of charge ordering phenomena. Different ground 
state phases may be realized depending on the strength 
on the Coulomb interaction. The details of these solu- 
tions will be described in the next subsections. 

Fig. |^ shows a representative phase diagram containing 
a cut along the V-S plane of the parameter space for fixed 
t/J = 1.25. Phase diagrams for other parameter values 
are shown in Figs. [6] and |?]. The inclusion of additional 
model parameters as a biquadratic exchange interaction 
can lead to furt her g round state phases, these will be 
described in Sec. 
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Let us start the discussion of Fig. |^ with large dop- 
ing. Here no phase separation tendency is present, the 
ground state charge distribution is homogeneous. The 
magnetic interaction together with the (infinite) on-site 
repulsion leads to pairing in the d x i_ y ? channel. Moving 
to smaller doping, the system encounters a phase sepa- 
ration instability at V = 0, i.e., the holes are expelled 
because the magnetic interaction favors a fully dimerized 
half-filled configuration. The Coulomb repulsion-counter- 
acts (frustrates) this phase separation tendency^ leading 
to microscopic charge ordering below a certain hole con- 
centration, S < <5stripc- The important point now is that 
the kinetic energy disfavors crystal-like states where the 
holes are essentially localized. Instead, the holes form 
"stripes" where hopping in one dimension is still pos- 
sible. Furthermore, the tendencies to fermion pairing 
on one hand and to dimerization on the other hand are 
still present. Consistent with this picture, our numeri- 
cal search always yielded lowest energy states with C and 
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FIG. 3. Ground state phase diagram of TL at N = oo, 
t/J = 1.25. Except for the d-wave superconductor all states 
have C broken. The "stripe" states have coexisting charge 
density wave order and superconductivity; the crystal states 
are insulating. The thick lines denote phase transitions being 
first-order in the large- A limit. Within the charge-ordered 
phases there are numerous additional transitions at which the 
detailed nature of the C symmetry breaking changes - these 
are not shown. The left and right boundaries of the stripe 
phase define 5pc(V) and 5stri pe (V), see text. 



5 broken, in a region 5pc < S < <5stripe and for small 
values of the Coulomb repulsion V. Th/ppe states consist 
of bond- centered charge density wavesL3'U (stripes) which 
co-exist with superconductivity. The effect of the kinetic 
energy becomes less pronounced when the doping is fur- 
ther decreased. Eventually, for extremely small doping, 

6 < Spc{t, J, V), and non-zero repulsion V the lowest en- 
ergy state is an insulating crystal of Cooper pairs which 
breaks C symmetry. Such a state arises from the com- 
bination of charge ordering and pairing tendencies. De- 
pending upon the parameter values t, t', J and V also a 
doped spin-Peierls state may be realized at small and/or 
intermediate doping. For strong Coulomb repulsion be- 
tween the charges (large V region in Fig. [^) one expects 
that hopping as well as pairing (mediated by the mag- 
netic interaction) will become unimportant. Then the 
situation resembles a low-density electron gas (where the 
potential energy dominates over the kinetic energy) , and 
the ground state becomes a Wigner crystal built of single 
charges. 

We continue with a more detailed description of the 
characteristics of the mean-field ground states. 



1. Homogeneous superconductor 



At doping levels smaller than a critical Sps(t,J) the 
"bare" t -j, J model in the large- A limit shows phase 
separation^ as explained above. [Note Sps(t, J) = 
^Stripe (t, J, V = 0).] In contrast, for 5 > 5ps and small 
V the ground states have uniform site charge distribu- 
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tions. These states include a doped spin-Peierls state 
(2x1 unit cell, see below) at relatively small 5, and 
states with a single site per unit cell which are homo- 
geneous superconductors. They can be characterized by 
two link field values Q x and Q y . In particular, one finds a 
d-wave superconductor with Q x = —Q y at intermediate 
5 (~ 20-50%), and an extended s*-wave superconductor 
with Q x — Q y for larger S. The excitation spectrum of 
the d-wave superconductor has four nodes whereas the 
s*-wave phase is fully gapped. 

Note that for t/J < 0.3 the critical doping level i5ps 
is in fact zero, so Titj does not show phase separation 
(and no stripes for V > 0) at very small t/J. Also, in 
the regime of small t/J the d-wave superconductor is re- 
placed by a state with s* + id x 2_ y 2 symmetry, i.e., the 
link fields obey Q y — Q x e l8 with a continuously vary- 
ing phase 9. Here the quasiparticle spectrum is again 
fully gapped, furthermore time-reversal symmetry T is 
broken. 

The described results for the homogeneous states are 
identical to the ones obtained in Ref. ||. Similar states are 
also found for large values of \t' /t\ where the tendency to 
phase separation is suppressed, see Sec. |IIB . 

2. Spin-Peierls state 

The spin-Peierls state which breaks C is found as 
ground state of the undoped system in the large- N limit. 
For 5 = its energy per site is given by E SP /(NN S ) = 
— J/4. Only one of the four link fields of the 2x1 unit cell 
is non-zero, i.e., the square lattice is completely covered 
by dimers. 

The corresponding doped state is the lowest-energy 
state with a homogeneous site-charge distribution in the 
small doping regimecl. The link fields in x direction 
Qix, Qix have different magnitudes reflecting the dimer- 
ization, whereas Q\ y — Q2y The fields cannot be made 
real simultaneously in any gauge, so the doped state 
breaks time-reversal invariance T and is fully gapped. 
However, if we restrict our attention to states which 
preserve T, then there exists a region of small dimer- 
ization, \Q\ X — Q2 X \ <C IQizI, near the transition to a 
d-wave state where the spectrum has 4 gapless points. 
Note that this state which has coexisting superconduct- 
ing and spin-Peierls order can be considered as a 2 x 1 
bond charge density wave with the hole densities being 
equal on all sit es_ Jlelated states have been studied in 
other approachesoEj, but they find spin-Peierls order in 
a Fermi liquid, not a superconductor. 

3. States with stripe charge order 

The saddle-point solutions at small doping, 5pc < 5 < 
^Stripe) and not too large V, break lattice translation sym- 
metry C and can be described as bond-centered charge 
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FIG. 4. The charge-ordering wavevector, K, (in reciprocal 
lattice units) as a function of S at N = oo for t/J = 1.25 
(Fig. []), and two different strengths of the Coulomb repulsion 
V/t = 0.5 and 1.0. For the stripe states as illustrated in Fig. [l] 
we have K = 1/p, the spin-Peierls state has K = 1/2 at 8 — 0. 
The K — 1 value at large 8 has C symmetry restored, and is 
a pure d-wave superconductor. For very small 8, the ground 
state is a Wigner crystal of Cooper pairs with incommensurate 
charge order. The unit cells for the largest plateaus are also 
shown. 

density waves. These states have a p x 1 unit cell, as 
shown in Fig. [l| We always found p to be an even in- 
teger, reflecting the dimerization tendency of the (5 = 
solution. 

The holes prefer to segregate in one-dimensional 
striped structures, i.e., within each p x 1 unit cell, the 
holes are concentrated on a q x 1 region. The link fields 
Q in the hole-rich region can be made real and have differ- 
ent signs in x and y direction reminiscent of <i-wave pair- 
ing correlation in the stripes. Most of the stripe phases 
are "fully" formed stripes, i.e., the regions between the 
stripes with a width of p — q have a hole density of zero 
and are insulating and fully dimerized. In this case the 
link fields between the stripe sites and the boundary sites 
of the insulating region vanish. However, for a parame- 
ter region close to the transition to a doped spin-Peierls 
state (see Fig. ^|) there occur "partially" formed stripes as 
ground states. These states have the same kind of charge 
density modulation as described above, but the hole den- 
sity in the hole-poor region is not zero, so the system is 
an anisotropic 2D superconductor even at iV = oo. In 
general, the hole density pi per unit length of each stripe 
is not of order unity as found in earlier theories but sig- 
nificantly smaller. The values of q and pi are determined 
primarily by t, J, V; they depend only weakly on S. For 
intermediate values of V where q = 2 we found values 







of pi ~ 1/2, here the stripe can be viewed as a ladder 
with roughly 1/4 hole per site. In general, smaller values 
of V yield larger values of q; the limit V — > leads to 
q — > oo which reflects the tendency to phase separation 
in the "bare" t — J model. 

The main effect of varying the total hole density is a 
change of p and therefore of the stripe distance. Wc.notc 
that the stripes do not survive in the limit 8 — > OEH: for 
8 < Spc the ground state changes to an insulating Wigner 
crystal of Cooper pairs. However, for small V, Spc is very 
small, <5 PC ~ exp(-l/V). For <5 PC < 8 < Stripe, we 
find an approximate proportionality p ~ l/<5, which also 
implies that pi is nearly independent of 8. The evolution 
of the ordering wavevector K = 1/p with 8 is shown 
in Figs. [|, ^, (J there are plateaus for each even integer 
number p. Our large- N theory only found "stripe" states 
in which K was quantized at the rational plateaus in 
Fig. ||. The reason for p, q being even and for the plateaus 
is easily identified as the strong dimerization tendency 
of the system. The columnar arrangement of the spin- 
Peierls singlets immediately leads to the "staircase" -like 
curve shown in Figs. f|, [|, ^. However, for smaller N we 
expect that irrational, incommensurate, values of K will 
appear, and interpolate smoothly between the plateau 
regions. 

For most values of t, J, and V there is a large plateau 
at p = 4 around doping 8 = 1/8, and, for some parameter 
regimes, this is the last state before C is restored at large 
(5; indeed p — 4 is the smallest value of p for which our 
mean-field theory has, solutions with hi not spatially uni- 
form. ExperimentallyHS, a pinning of the charge order 
at a wavevector K — 1/4 is observed, and we consider 
it significant that this value emerges naturally from our 
theory. 

It is worth pointing out that the hole density pg per 
unit length of stripe is not exactly pinned at one value 
which means that the stripes are not incompressible. The 
hole density varies continuously within each plateau but 
jumps discontinuously as the transition is made from one 
plateau to the next. 

There are strong pairing correlations between the holes 
in each g-width region. Strictly speaking, for N = oo 
each q- width stripe above is a one-dimensional supercon- 
ductor, while the intervening (q — p)-width regions are 
insulating. However, fluctuation corrections will couple 
with superconducting regions, yielding an effective theory 
discussed in Section VII of Ref. [54] with their dimension- 
less parameter K ~ N. This implies that Josephson 
pair tunneling between the one-dimensional supercon- 
ductors is a relevant perturbation at sufficiently large N, 
leading to a two-dimensional anisotropic superconduct- 
ing ground state, and should allow good metallic conduc- 
tion above the superconducting transition temperature^ 
These characteristics are consistent with observationsEj 
on La2-x-j/Ndj,Sr x Cu04. 

The quasiparticle spectrum in the striped phases is al- 
ways fully gapped because of the presence of supercon- 
ductivity. Except for the partially striped phases where 
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FIG. 5. Fermionic mean-field excitation spectrum in fully 
striped phases, t/J = 1.25. The spectrum is independent of 
k x , the momentum perpendicular to the stripes. Left: q — 2 
phase with p# = 0.8. Right: q = 4 phase with pi = 0.6. In 
the second case the bandwidth is much smaller due to the 
smaller doping level in the stripes (pi/q). The flat bands at 
energy J/2 correspond to excitations of localized dimers in 
the undoped regions. 

the fermion spectrum shows a weak dependence on the 
momentum k x perpendicular to the stripesQ, the spec- 
trum of the stripe phases is fc^-independent. The fermion 
spectrum in the partially striped phases were displayed 
in Ref. ^, and in Fig. ^| we show the quasiparticle bands 
for the fully striped cases. The spectrum consists of two 
contributions: The excitations of the fully dimerized un- 
doped regions correspond to removing one fermion from 
a localized dimer, the energy cost is J/2 and of course 
momentum-independent. The stripes of width q give rise 
to q one-dimensional dispersing bands where the gap is 
determined by the pairing amplitude of order J, and the 
total bandwidth is given by Atp?/q where pi is the linear 
hole density per unit stripe as above. The stripe exci- 
tations therefore lie in the gap of the spin-Peierls state; 
their properties are entirely determined by pi and q (and 
the model parameters t and J). Upon increasing the 
doping level 8 within one of the plateaus in Fig. [I] the su- 
perconducting gap decreases slightly since the link fields 
Q decrease, and the position of the dispersion minimum 
changes corresponding to the Fermi momentum. The 
transitions between the plateaus arise from level cross- 
ings, i.e., there is no critical dynamics associated with 
them. At the transition from one plateau to the next, pi 
decreases discontinuously, leading to an increase in the 
gap and so on. If the stripe width q changes at a plateau 
transition then the number of dispersing bands will also 
change, see Fig. ||. 

At this p^ip^ we briefly discuss earlier mean-field 
calculations£3c3 which predicted an inhomgeneous 
charge distribution in the ground state of the Hubbard 
and related models. These computations where based on 
the observation that large-,!? mean-field theories of doped 
antiferromagnets (S denotes the size of the spin) show 
charge density wave instabilities - the same applies to 
the large-iV theory discussed here. An important dif- 
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FIG. 6. Same as in Figs. | and §, but for t/J = 0.5. 
Upper panel: Ground state phases of 7i as function of dop- 
ing and Coulomb repulsion. Thick and thin lines denote 
first-order and second-order transitions, respectively. Lower 
panel: Charge-ordering wavevector, K, (in reciprocal lattice 
units) for V/t = 0.2. For 28% < S < 42% the ground state 
is the doped spin-Peierls state breaking C and S; it has a 
uniform site-charge distribution and a 2 x 1 unit cell. 



ference, however, is the character of the stripes. r-Ttci 
early mean-field calculations for the Hubbard modecHTLJ 
predict insulating, site-centered stripes with a hole den- 
sity of unity within the hole-rich regions; the experimen- 
tally found stripe states in the cuprates are, however, 
either metallic or superconducting. In contrast to these 
early mean- field results, the present large- N computa- 
tions yield superconducting, bond-centered stripes, and 
it would be useful for future experiments to detect this 
distinction between bond- and site-centering. Let us dis- 
cuss this difference more precisely for the case p = 4. For 
the site-centered stripe, the hole density per unit length 
in each column of sites takes values pi, p2, P3, pi be- 
fore repeating periodically (where pi_3 are some three 
distinct densities); this is the configuration usually as- 
sumed in most experimental papers. In contrast, for the 
bond-centered state, these densities take the values pi, 
Pi i p2i P2, and this also appears to be compatible with 
existing observations. The bond-centering is important 
in the present computation for enhancing pairing correla- 
tions, which are responsible for superconducting/metallic 
transport in the direction parallel to the stripes. 

Possible differences in the spin dynamics between 
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FIG. 7. Same as in Figs. | and % but for t/J = 2.5. 
Upper panel: Ground state phases of Ti as function of dop- 
ing and Coulomb repulsion. Lower panel: Charge-ordering 
wavevector, K, (in reciprocal lattice units) for V/t = 1.0. 



bond- and site-centered striped phases have been dis- 
cussed by Tworzydlo et alu. They argued that the mag- 
netic domains between the stripes can be described by 
spin ladders with either an even or odd number of legs; 
these two cases lead to very different spin fluctuation 
spectra. _ , 

We note that a very recent NQR experimenf£3 indi- 
cates a charge distribution which has some features con- 
sistent with the bond-centered state: they find only two 
inequivalent sites associated with the stripes, and a den- 
sity in the hole-rich region which is considerably smaller 
than 0.5. 



4- Wigner crystals of Cooper pairs 

For very small doping the "stripe" states become un- 
stable with respect to further C breaking. This leads 
to two-dimensional insulating states with S restored, 
whiclt-can be characterized as Wigner crystals of Cooper 
pairsEd. In the large- N limit, the Cooper pair crystals 
consist of pairs of adjacent sites with non-zero hole den- 
sity and a non-zero bond amplitude Q between them; the 
remaining sites of the square lattice form singlet bonds as 
in the fully dimerized, undoped spin-Peierls state. The 
occurrence of such pair crystal states is intimately re- 
lated to the dimerization tendency of the system: In a 
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pair crystal the repulsion energy is much smaller than a 
stripe state, and the hole-rich sites still form dimers lead- 
ing to a non-zero contribution from TLtj- (At N = oo the 
hole density on the sites of each pair is not exactly unity, 
b\ < 1, this is an artifact of the large- AT limit. For phys- 
ical values of N the physical system will certainly have 
two electrons per pair.) 

To minimize the repulsion energy, the hole pairs want 
to form a triangular Wigner lattice (with a lattice con- 
stant ~ V~5). However, in the present problem the hole 
positions have to be on the underlying square lattice. 
This leads to incommensurate structures where the pairs 
form an approximate triangular lattice. The energy per 
site in this state at small 6 is given by Epc/(NN S ) = 
—J/4 — a 6 where a contains contributions from t, J, 
and V. 

We note here that there are also crystal-like insulating 
saddle-point solutions with clusters of sites with non-zero 
hole density. These "charge islands" can be considered 
as an integer number of Cooper pairs bound together - 
the described pair crystals are the smallest possible is- 
lands. Larger islands occur at small doping, they replace 
the pair crystal phases shown in the phase diagrams at 
small hopping t. We have not systematically examined 
all possible "island" phases because large islands require 
large unit cells in the numerical calculation. However, 
by estimating the ground state energy of the "island" 
phases we have checked that stripe phases are always fa- 
vored in the experimentally interesting doping regime, 
8% < S < 20%. 

Summarizing these findings, the frustration of com- 
plete phase separation by V leads to two possible sce- 
narios of charge clusters with "smaller" size: (i) stripes 
which are Id objects, i.e., the charges are confined in one 
direction (note that finite T may melt the stripe order 
leaving behind finite length parts of stripes; strong im- 
purity influence may also break up stripes into segments), 
and (ii) islands which can be considered as bound states 
of Cooper pairs - here the charges are confined in both 
directions. The calculations show that for physical values 
of t/J and not too small doping S the stripe scenario (i) 
is favored due to the larger kinetic energy of the holes in 
the stripes. 

5. Wigner crystals of single charges 

For large V the Coulomb repulsion dominates over t 
and J which leads to the appearance of Wigner crystal 
states - here the crystal consists of single charges. The 
holes minimize the Coulomb energy by forming a crystal- 
like structure - this is similar to the charge ordering in 
a low-density electron gas. The lowest Coulomb energy 
is again reached for the case where the charges form a 
triangular lattice with a Wigner lattice constant ~ \J~8~, 
however, the hole positions have to coincide with the sites 
of the underlying square lattice. Therefore, the states 



at large Coulomb repulsion V will have structures in- 
commensurate with the square lattice, i.e., a (infinitely) 
large unit cell with holes sitting on a fraction 1/6 of the 
sites (b 2 = 1). These sites form an approximate triangu- 
lar Wigner lattice, all other sites are half-filled (b 2 = 0) 
and form dimerized bonds. (The calculation shows that 
at finite values of V/ J the competition of exchange and 
Coulomb energies at TV = oo leads to hole densities on 
the sites of the Wigner lattice being smaller than unity, 
b 2 < 1, and to a distortion of the ideal dimer pattern in 
the link fields Q near each hole site. This again is an 
artifact of the condensation of the slave bosons in the 
large- A limit.) 

The kinetic energy in the Wigner crystal state van- 
ishes. The exchange energy per site is given by — J/4 x 
(1 — 5) at low 6 since (1 — 6) is the number of sites remain- 
ing for singlet formation. It is important to note that the 
Coulomb energy is negative, i.e., {TLv) ~ — ^ S 1 ^ 2 < 
with 7 ~ V. So we can estimate the energy per site for 
small 6 to be E WG /(NN S ) = -J/4(l - 6) - 7 <5 3 / 2 . 

B. Influence of additional model parameters 

The calculations presented so far show that the ten- 
dency to stripe formation at low doping is a very robust 
feature of the t — J — V model at N = oo. In the fol- 
lowing we discuss the influence of the inclusion of further 
processes into the Hamiltonian. 

1. Longer-range hopping 

We start with an additional next-nearest neighbor hop- 
ping t' which is supposed to be important in certain 
cuprate superconductors. The large- A calculations show 
that the results are nearly independent of t' as long as 
it is small, \t' /t\ < 0.3. Larger values of t' of either sign 
suppress the tendency to phase separation in the "bare" 
model (without TLv)- This in turn means that the stripe 
instability is weakened; the stripe region in the phase 
diagram shrinks, it is shifted to lower doping and lower 
values of V. This is .consistent with e.g. the findings of 
White and ScalapincEa. 

2. Biquadratic exchange 

We have also considered the addition of a biquadratic 
spin exchange 

H * = §s E^-S,) 2 (2.9) 

<hj> 

to the Hamiltonian fl2.4|) for general N. The strength of 
the biquadratic interaction is given by k J. For the phys- 
ical case N = 1 the biquadratic exchange term simply 
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reduces to S, • Sj + const, but for larger A^j-jJb contains 
exchange processes involving eight fermionsEj. For the 
total coupling of two (isolated) spins to be antiferromag- 
netic there^ an upper limit for k which is given by k < 2 
for any Alia. In the large- N limit the decoupling can be 
done in two steps, it turns out that no new link fields are 
necessary. At the saddle point the additional term takes 
the form: 



Ha = 



H.J 



E 

<i,j> 



\Q 4j \ 2 (2J^flf} p Q* 3 + H.c. - 3N\Q l3 \'- 



(2.10) 

As discussed in Ref. |l^ such a biquadratic exchange term 
weakens the tendency to dimerization. In fact, in the 
SU(iV) saddle-point approach it can be used to stabilize 
a staggered flux phase; it is expected to find a similar 
effect on the hole-poor regions of the stripe statesEII. 

In the undoped system, we find that the N = oo 
ground state changes from spin-Peierls (at k = 0) to 
a box phase for < k < 1.13, and then to a "flux"- 
like phase for k > 1.45. (Some intermediate phase oc- 
curs for 1.13 < k < 1.45.) All phases at J > can be 
described within a 2 x 2 unit cell. The box phase has 
non-zero link fields Q around isolated-pJaquettes of four 
sites each, all other link fields vanishEj. In the "flux"- 
like phase all Q are equal in magnitude, their phases can 
be made real by a gauge transformation, and the expec- 
tation value of the plaquette operator IIq is real and 
negative. However, unlike in the SU(7V) case there is 
no direction associated with the Q link fields. A direct 
interpretation in terms of a flux is therefore not possi- 
ble. The described mean-field solutions are identical with 
the results p£ Marston and Affleck for the SU(7V) Heisen- 
berg modelEj. In fact, by applying a generalized particle- 
hole transformation it can be shown that the mean-field 
equations for the undoped system are equivalent for both 
SU(2A0 and Sp(27V) cases. 

Turning to finite doping, we find that all homogeneous 
phases are unstable w.r.t. phase separation at V = 
and low doping. The inclusion of V l eads to stripe states 
similar to the ones described in Sec. II A 3 for all values 
of k. The hole concentration pi in the stripes is again 
determined by the microscopic parameters, and is found 
to be between 0.2 and 0.5 for reasonable values of t, J, 
and V. The stripes are one-dimensional superconduc- 
tors, for most of the parameter values the bond ampli- 
tudes Q within the stripes can be made real and reflect 
d-wave pairing correlations [at larger V and intermediate 
6 there is a small region with (s* + id)-likc pairing in the 
stripes]. The Q fields in the insulating regions between 
the stripes reflect the structure of the undoped ground 
state: At small non-zero k one finds a box-like structure 
whereas k > 1.45 gives rise to a more homogeneous dis- 
tribution of the Q which can be understood as flux-like 
arrangement with distortions at the boundaries (stripes). 
Interestingly, in the latter case there can occur non-zero 
link fields Q between the stripes and the undoped region 
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FIG. 8. Ground state phase diagram for the Sp(2iV) t — J 
model with biquadratic exchange term as function of doping 
and Coulomb repulsion, t/J = 0.75, k = 1.5. 



(leading to longer-ranged spin correlations), but these 
disappear for larger hole concentration in the stripes. 

At larger doping the stripes are disfavored, the ground 
state becomes a homogeneous superconductor with either 
d-wave or s* +id symmetry. Compared to the case k = 0, 
the tendency to phase separation is slightly weakened by 
the introduction of the biquadratic exchange, the region 
of stripe ground states shrinks. A sample phase diagram 
is shown in Fig. ||. Note that the transition from d-wave 
to s* + id symmetry with decreasing doping is accompa- 
nied with the opening of a gap in the fermionic spectrum 
(and with spontaneous T breaking). 



3. Next-nearest neighbor exchange 

To explore further possible symmetries of the super- 
conducting order parameter (which is represented by 
the link fields Q in the mean-field theory) we con- 
sider an additional next-nearest neighbor exchange J' > 
0. Note that the inclusion of such a frustrating next- 
nearest neighbor exchange between pairs of sites on the 
same sublattice is straightforward within the present 
Sp(2./V) mean-field theory since the same representation 
of Sp(27V) is employed for spins on both sublattices [in 
contrast to the bosonic SU(N) approach]. The new in- 
teraction is decoupled similar to the nearest-neighbor in- 
teraction as described at the beginning of Sec. ||; this 
leads to diagonal link fields Q\y The model with hop- 
ping terms t and t' , exchange processes J and J' as well as 
biquadratic exchange strengths k J and k J' shows a num- 
ber of new ground state phases at the mean-field level, 
among them superconductors with mixed order parame- 
ters, spin-Peierls and flux-like phases as well as charge- 
ordered states. 

As an example, we show the phase diagram for t/ J = 
1.0, J' I J = 0.4, and k = 1.5 in Fig. §. The un- 
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FIG. 9. Ground state phase diagram for the Sp(2iV) 
t — J — J' — V model with biquadratic exchange term as func- 
tion of doping and Coulomb repulsion; t/J = 1.0, J'/ J = 0.4, 
t' = 0, K = 1.5. d + id denotes the T-breaking superconductor 
with (d x 2_ y 2 + id xy )-wecve symmetry. 



doped system has a large ground state degeneracy, small 
doping produces stripe states due to the phase separa- 
tion instability in the absence of Couplomb repulsion 
V. The stripes disappear at larger doping, leading to 
a superconductor with a homogeneous charge distribu- 
tion. First we find a (d x 2_ y 2 + id xy )-w&ve superconduc- 
tor; further increasing 5 leads to a pure d x 2_ y 2 state. The 
{d x 2_ y 2 + id xy )-w&ve state has a fully gapped spectrum 
and breaks time reversal symmetry. It is characterized by 
four Q link fields which obey Q x = —Q y , Q'n = —Q[ _i 
and Q'n = ieQ x . Here, e is a real value measuring the 
mixing of both d-wave order parameters. At the transi- 
tion line between the (d x 2_ y 2 + id xy ) state and the pure 



-wave state the e vanishes continuously. 



4- Alternative interaction decoupling 

There has been recent interestSil in possible ground 
states which spontaneously break timepceversal symme- 
try T and show circulating currentscjii3 (doped "flux" 
phases) . Such states can be characterized as orbital mag- 
nets; if the current direction around each unit cell alter- 
nates in sign between adjacent cells one obtains an orbital 
antiferromagnet. To explore possible realizations of such 
states within our saddle-point approach we extend the 
Hamiltonian by an interaction term which provides de- 
coupling in the particle- hole channel in the large- N limit. 
It has the form (- J tf )/(2JV) {flfjaXfLfif}) which also 
reduces to Sj • Sj + const for N = 1. The interaction is 
decoupled by link fields x which take the values 

NXij = Nx% = (fLf ja ) (2.11) 

at the saddle point. In contrast to the Q fields these new 
fields Xij specify a direction for the link (ij). For the 



undoped case the phases of the x are gauge-dependent, 
but the plaquette operator II X = X12X23X34X41 is again 
a gauge-invariant object. In the doped system, the fields 
X are directly connected with the current transported 
through one link, jij = 2NIm.Xij, so their phases are 
meaningful. 

At the saddle point the Hamiltonian containing the 
interaction J takes the form 



i>j 



J, 



//a/^+H.C.-iV|^| : 



(2.12) 



which is formally a contribution to the fermion hopping. 
We also include a corresponding biquadratic interaction 
which can be decoupled in the particle-hole channel lead- 
ing to 



Hi 



H.J 



<i,j> 



| 2 (2/i/ JQ X*,+H.c.-3iV|x y | 



(2.13) 



The mean-field eq uatio ns show that a magnetic interac- 
tion of the type ( 2.12j ) is equivalent to the decoupling 
done within a SU(2A) large-iV approach to the anti- 
ferromagnet. The sum of both terms, Tij + Hj, con- 
tains both possible factorizations (particle-particle and 
particle-hole) of the four-fermion interaction term. It 
is therefore equivalent to an unrestricted Hartree-Fock 
treatment of the original Heisenberg Hamiltonian for 
N = 1 written in terms of auxiliary fermions, provided 
that Jij = Jij . In the following, we will discuss the large- 
TV limit and treat J and J as independent parameters, 
concentrating on the region where J ~ J. Varying the 
ratio x = J/(J+ J) can be viewed as continuous interpo- 
lation between the Sp(2A) (x = 0) and SU(2iV) (x = 1) 
approaches to the antiferromagnet. 

For x = 1/2 (J = J) and doping 5 — 0, we find a 
large ground state degeneracy which has also been re- 
ported earlier in unrestricted Hp;tree-Fock treatment of 
the Heisenberg antiferromagnetEH. The degeneracy can 
be lifted by the introduction of a biquadratic exchange 
k > which leads to box-like phases having a 2 x 2 unit 
cell. These states have non-zero values of both the Q and 
X link fields. 

Turning to finite doping, we find no phase separation 
instability at x = 1/2 independent of k — at zero or small 
V states with a homogeneous site charge distribution al- 
ways have lower energy than inhomogeneous states; there 
are no stripe states at x = 1/2. The ground state at low 
doping is either a spin-Peierls state (at k — 0) or box- 
like (k > 0), for larger doping it becomes a pure d-wave 
superconductor. The phase diagram for k — is shown 
in Fig. |l0|. The spin-Peierls state shows dimerization in 
both th e Q an d the x fields with a pattern as described 
in Sec. 



IIA2 



It is interesting to note that the doped 
spin-Peierls state does not break time reversal symmetry 
T, so the Q fields can be made real by a gauge trans- 
formation (in contrast to the spin-Peierls state occurring 
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FIG. 10. Ground state phase diagram for the t — J model 
with interactions Tij and Hj, t/J = 1, x = 1/2 (J = J), 
k = (no biquadratic exchange). There are two spin-Peierls 
phases shown: At larger doping the dimerization is small lead- 
ing to gapless spectrum with four nodes. Upon decreasing 5 
the dimerization increases, and when it exceeds some critical 
value, the spectrum becomes gapped: the gap opens when the 
nodal points collide in pairs at points on the lines k y = ±7r/2 
in the extended Brillouin zone. 
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FIG. 11. Ground state phase diagram for the t — J model 
with interactions Tij and Tij which are both supplemented 
by a biquadratic exchange, k = 1.5. The remaining param- 
eters are t/J = 1.25, x = 0.6 (J = 1.5J). OAF denotes the 
T-breaking orbital antiferromagnet with staggered circulat- 
ing currents around each plaquette. At S — we have the 
staggered flux state. 

in the "pure" Sp(2/V) case which is always T breaking). 
Similarly, the doped box state for x = 1/2 also preserves 
T. 

As one further example, we describe the ground states 
at x = 0.6. The undoped system without biquadratic 
exchange (k = 0) again shows a large degeneracy. Values 
of < k < 1.4 lead to a box- like phase. For k > 1.4 we 
find the staggered flux phase as ground state, here all \ 
are equal in magnitude and the plaquette operators IT X 



have the negative real expectation values. The particle- 
particle link fields Q vanish for any k. 

At finite doping, there is again no phase separation 
tendency at x = 0.6. For k = the ground state phases 
are spin-Peierls at small S and d-wave at larger S, whereas 
small non-zero k produces box-like ground states. Inter- 
esting phases occur for n > 1.4 as shown in the phase 
diagram in Fig. [ll| The staggered flux phase evolves 
to an orbital antiferromagnet with a two-site unit cell 
and circulating currents (given by Imv) in each plaque- 
tte with alternating direction betwees, adjacent plaque- 
ttes. This state is unstable to pairingt-3: the Q fields ac- 
quire non-zero, real values leading to a T-breaking state 
with coexisting orbital antiferromagnetism and d-wave 
superconductivity!^. We note that this coexistence oc- 
curs here at the mean- field .level, this can be contrasted 
to other mean-field theoriesO where the flux and d-wave 
instabilities occur in different orders in a 1/N expansion. 

For larger doping, the currents vanish, and the ground 
state becomes a pure d-wave superconductor. All the de- 
scribed states have a gapless excitation spectrum - note 
here that the flux phase can be interpreted as a d x 2_ y 2 
density wave, implying nodes of the order parameter in 
the diagonal directions in momentum space. 

At very large values of V the dominating repulsion 
gives rise to Wigner crystal states as described above. 
However, there are no pair crystal phases for the cases 
x = 0.5 and 0.6 which is related to the absence of phase 
separation tendencies at these parameter values. 

Finally we note that the case x — 1, i.e., the "pure" 
SU(iV) interaction decoupling, leads to pitase separation 
in the absence of Coulomb repulsion VEj, and the in- 
clusioiLof V is expected to produce "striped" ground 
statesO. 



C. Phase transitions at TV = oo 

We briefly discuss the nature of the transitions between 
the various ground state phases at the mean-field level 
(N = oo). 

The transition from a d-wave superconductor, with 
C unbroken, to the fully-formed p x 1 stripes discussed 
above, can either be first-order, or via intermediate states 
with partial stripe order. In the latter case, there is first a 
continuous transition to a doped spin-Peierls state with C 
symmetry breaking at p — 2. To our knowledge a p = 2 
charge-ordered superconducting state has not been ex- 
perimentally detected, but a search for one should be 
worthwhile. Then there is a second second-order transi- 
tion to p — 4 state with partial stripe order, before the 
fully-formed p = 4, q = 2 state with intervening insu- 
lating stripes appears. In any case, the excitation spec- 
trum in the stripe states is fully gapped. The transitions 
between the different stripe phases where p changes its 
value are in general first-order since they arise from level 
crossings. 
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At very small hole densities the Cooper pair crys- 
tal states are favored over striped states, then the 
ground state nature changes from one-dimensional to 
two-dimensional C breaking. Again, this transition is ei- 
ther first-order or second-order depending on the values 
of t, J, and V, in the latter case it can be visualized as 
additional C breaking in the longitudinal stripe direction. 

Finally, the transition to a Wigner crystal phase of sin- 
gle charges at large V is always a first-order transition 
because the inherent tendency to fermion pairing is bro- 
ken in the Wigner crystal. 

By adding a biquadratic exchange term to the Hamil- 
tonian we found a T-breaking superconductor with s*+id 
symmetry and a gapped spectrum. As displayed in the 
phase diagram of Fig. ||, it shows a second-order transi- 
tion to a pure d-wave state with increasing doping. Sim- 
ilarly, the introduction of a second-nearest neighbor ex- 
change J' can lead to a T-breaking (d x 2_ y 2 + id X y)-wave 
superconductor which also has a second-order transition 
to a pure d-wave state with increasing 8, see Fig. O. 

The SU(2A0-like interaction term Hj (Sec. pB4|) 



A. d-wave to (s* + id)-wave or 



+ id x 



which has to be decoupled in the particle-hole channel 
leads to other possible transition scenarios. As shown 
in Fig. |l^, a T-conserving spin-Peierls phase is possible. 
Then two transitions are found upon decreasing doping: 
First there is a second-order transition from the d-wave 
state to a (weakly dimerized) spin-Peierls state where 
the dimerization sets in. The four nodes of the d-wave 
phase are preserved, but shift gradually in momentum 
space. At smaller 8 and a "critical" dimerization there 
is another second-order transition where the spectrum 
becomes gapped. Turning to the x = 0.6 phase dia- 
gram of Fig. [ll], we note that the transition between 
the OAF+d-wave coexistence phase and the pure d-wave 
state is second-order. Both order parameters have d x i_ y i 
symmetry leading to gapless spectra. The T symmetry 
is broken in the low doping phase with orbital currents. 



III. CONTINUOUS QUANTUM PHASE 
TRANSITIONS 



In Section II C we noted a number of continuous quan- 



tum phase transitions in the mean-field phase diagrams 
of Fig. p[-pd|. As we discussed in Section 1 A , for the pur- 



poses of understanding the photoemission experiments 
of Ref M we are particularly interested in transitions 
in class A. Of these, the simplest involve T symmetry 
breaking alone, and we will consider these in the first sub- 
section below. Transitions involving breakin g of C sym- 
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This 



metry alone are considered next in Section 
is followed by a discussion of the onset of staggered flux 
order in a d-wave superconducto r (thi s transition breaks 
T and C symmetries) in Section 1IIC. 



A continuous transition between a d-wave supercon- 
ductor and a (s* + id)-wave superconductor appears in 
the mean-field phase diagram in Fig. and fluctuations 
in its vicinity will be described below. A very closely re- 
lated theory applies to the transition between a d-wave 
superconductor and a (d x i_ y i + id xy )-wave superconduc- 
tor shown in Fig M. Indeed, to the order we are comput- 
ing things, the results for the scaling functions of the two 
transitions are identical. 

We note that a transition between d and (d x 2_ y 2-L£d xy ) 
pairing was also considered recently by D.-H. LeeO us- 
ing a rather different framework. Lee concluded that the 
transition was in the universality class of the Ising model 
in a transverse field; this is equivalent t o as suming that 
the critical theory is described by S$ in ( |3.3| ) below, and 
places the transition in class B. In fact, as we show be- 
low, it is also necessary to inc lude fermionic excitations: 
the complete theory is in (3/7) below, and the transition 
is in a new universality class which belopgs in class A. 
In a separate recent work, Balatsky et cd.cil considered a 
"clapping" collective mode in a (d x 2_ y 2 +id xy ) supercon- 
ductor: this is an 'angular' fluctuation mode expected in 
a state with well-formed d xy pairing induced by an exter- 
nal magnetic field; in contrast we are interested here in 
spontaneous fluctuations to (d x 2_ y 2 + id xy ) order, which 
will dominated by 'amplitude' fluctuations represented 
by our field <f> below. 

We begin by establishing notation, and reviewing the 
familiar physics of the low-energy fermionic quasiparti- 
cles in a d-wave superconductor. For thc-jnost part, we 
will follow the notation of Balents et al£3. The low en- 
ergy excitations appear at four points in the Brillouin 
zone: (K,K), (-K, K), (-K,-K), and (K,-K). We 
denote the components of the electron annihilation oper- 
ator, c a , in their vicinity by f la , ha, ha, ha respectively, 
where a =|, J, is the electron spin component. Now intro- 
duce the 4-component Nambu spinors = (fi a , £abfl ) 
and ^2 = (ha, £abfl ) where E ao is an antisymmetric ten- 
sor with e^i = 1. The action for the fermionic excitations 
in the d-wave superconductor is then 



S-qr 



d d k 



(27T) 



-jT^2 *i (-««>„ + v F k x r z + v A k y r x ) #i 



A^jTV {-iu) n + v F k y r z + v A k x r x ) * 2 . (3.1) 



Here t" are Pauli matrices which act in the fermionic 
particle-hole space, k x>y measure the wavevector from the 
nodal points and have been rotated by 45 degrees from 
the axes of the square lattice, and vf, va are velocities. 

Let us now consider the transition to the (s* + id)- 
wave superconductor. This involves only a Z 2 symmetry 
breaking of the T symmetry, and so the order parameter 
is a real scalar field, <j). In the presence of a non-zero, 
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spacetime-independent (j>, the superconducting gap 
tion takes the form 

(c k -[C-ki) = A (cosfc 2: - cosky) + i 



func- 



cos k x + cos k y ). 



(3.2) 



The superconducting order parameter also has a single 
overall complex phase, but because the superfluid stiff- 
ness is finite, its fluctuations can be neglected in the crit- 
ical theory. On general symmetry grounds, we can write 
down the following effective action for <f> fluctuations 



S A = d d xdr 



(3.3) 



Here c is a velocity, sq is the parameter which tunes the 
system across the transition, and uq measures the quartic 
self-interaction of the order parameter. Normally in the 
absence of strict particle-hole symmetry, a term with a 
first-order time derivative can potentially appear in ; 
however, </> is a real field, and the only possible relevant 
term, (f>d r (f> = (l/2)<9 r </) 2 , is a total derivative and inte- 
grates to zero. 

The final piece of the action is the coupling between 
4> and the fermionic excitations. This can be easily de- 



duced from (3^), and the standard pairing interaction 
between A(fc) and the electrons. In the vicinity of the 
nodal points, this coupling reduces to the very simple 
expression^ (we are assximing here that K ^ n/2, as is 
the case experimentally^, so that (cosk x + cosfc y ) does 
not vanish at the nodal points) 



<S<i>0 = / d d xdr 



A, 



(3.4) 



The remainder of this subsection will consider the prop- 
erties of the theory + + 5$^. Before we embark 
on this case, we note the generalization to the case of 
a transition to d x 2_ y 2 + id xy order: in this case (3^2) is 
replaced by 

(cfcjc_fc|) = Aq(cos& x — cosky) + i(f>sink x sin k y ; (3.5) 



in the_sacinity of the nodal points, instead of (3.4), we 
obtainEl 



d d xdr 



(3.6) 



The only change is the relative minus sign between the 
two terms, which ari ses from the changing sign of the 
factor sinfca; sink y in (3.5) between the nodal points; it is 
this change which is responsible for thejipn-zero spin Hall 
conductance of the d x 2_ y 2 + id xy stateE3. The properties 
of the theory + S$ + are very closely related to 
that of Sqi + Scf, + S indeed, at the one-loop order we 
compute things here, there is no difference between the 



theories, and our results can be applied equally to the 
transition to cither s* + id or d x 2_ y 2 + id xy order. 

Our rcnormalization group analysis is closely related 
to that discussed by Zinn-JustirO. We will mainly re- 
strict our study to the case in which all the velocities 
in Sq,^ are equal to each other. By a suitable choice 
of the scale of time we can then set c = vp = v/\ = 1 
at the outset. Deviations from exactly equal velocities 
will be considered in Section HI A 1 below, and in future 
work. With exactly equal velocities, the action is actu- 
ally Lorentz invariant and so we must have the dynamic 
exponent z = 1; this will value of z will be implicitly as- 
sumed in the remainder of this paper. We introduce no- 
tation to make the Lorentz invariance explicit. We map 
*2 {t x +r z )* 2 /\/2, and define W x , 2 = -ir v ^\ t2 , and 
introduce Lorentz index /i = r,x,y. Then + + S^^ 
can be written in the compact Lorentz-invariant form 



S. 



s* -\-id 



= d 



(3.7) 



where D = d+ 1, and 7 P = (—t v ,t x ,t z ). We note that 
there is also a corresponding action Sd 



which 



(using (3.6)) differs only in having a relative plus sign 
between the couplings of 'J'i'I'i and ^2^2 to (j>. 

We now proceed with the renormalization group analy- 
sis of S s * + id or Sd 2 _ 2+idxy- We will present the results 
using the field-theoretic method as it leads to a more 
streamlined discussion of finite temperature properties. 
To this end, we introduce the wavefunction renormaliza- 
tions 



cj) = zl /2 4> R 

#1,2 = ^) /2 *l,2fl 

and the coupling constant renormalizations 



(3.8) 



An = 



u 



fi 



/2 



Zx 



jl/2 



z f z b 



1/2 



A 



J^_Zn 
SnZl 



(3.9) 



where fi is a renormalization momentum scale and Sd — 

One loop calculations of the renormalization constants 
using minimal subtraction of poles in e = 4 — D = 3 — d 
yield 



Z h = 



Z f = 



Z, 



Z u — 



4A 2 
e 

1-* 

2e 
A 2 

1 + — 

e 

3m 
1 + Ye~ 



48A 4 



(3.10) 
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From these we obtain the beta-functions (note that 
we are Hsing the field-theorists' definition of the beta 
functions, and these are opposite in sign to those nor- 
mally used by condensed matter physicists): 



/3(A) = -|A + ^A 3 
3 

(3(u) = -eu + -u 2 + 8u\ 2 - 48A 4 . 



(3.11) 



These equations have the infrared stable fixed point 



,*2 



7 

16e 
~21 



(3.12) 



which controls the physics in the vicinity of the quan- 
tum critical point. In particular, the fields acquire the 
anomalous dimensions 



Vb = 4A* 2 
rit = A* 2 /2 



(3.13) 



which will play an important role in the spectral func- 
tions to be discussed below. 

In the first subsection below, we will consider the con- 
sequences of unequal velocities. The next subsection will 
describe the T > quantum critical spectral response 
functions of the fixed point (3.13). 



1. Unequal velocities 

The full treatment of the problem with unequal veloc- 
ities is considerably more complicated than the simple 
renormalization group analysis above, and is deferred to 
future work. It is entirely possible that there are other 
non-Lorentz-invariant fixed points describing the criti- 
cal theory: however, apart from differences in the de- 
tails of the scaling functions, the qualitative properties of 
such fixed points should be quite similar to the Lorcntz - 
invariant case explored in more detail in Section III A 2. 

We will restrict our attention here to perturbations 
which break Lorentz invariance in a linear stability anal- 
ysis of the Lorentz- invariant fixed point. Such an analysis 
was also carried out in Ref. but the authors do not ap- 
pear to have performed the proper decomposition of the 
perturbations into the appropriate eigenoperators, as is 
required to obtain the correct scaling dimensions. 

It is important to classify the perturbations in terms 
of irreducible representations of the Lorentz symmetry 
of the fixed point. For the velocity differences, these are 
the symmetric, traceless, second-rank Lorentz tensors. 
So we consider the following perturbation to the action 
5<p + S<f, + SV^: 



S a 



d°X 



(3.14) 



It is crucial that the tensors above are traceless: this 
ensures, at linear order, that such perturbations do not 
mix with terms already in the Lorentz-invariant action, 
and with redundant operators which can be removed by 
a rescaling of the fields which are being integrated over. 
It is a simple matter to compute the one-loop renor- 



malization group flow equations to the terms in (3.14) 
We find 



/3(^i } ) 

0(9$) 



= A 



4 -2 -2 
1 1 -1/3 1/3 
1/3 1/3 




(3.15) 



Clearly, we need the eigen value s and eigenvectors of the 
matrix of coefficients in (3.15). First, note that there 
is an eigenvector, (1,1,1), with zero eigenvalue. This 
is expected, and its presence is an important check on 
our calculation — it corresponds to the perturbation of 
simply changing all velocities by the same amount: such 
a redundant perturbation can be absorbed by a rescaling 
of spacetime co-ordinates, and does not modify any of 
the physics. 

The other eigenvalues and eigenvectors correspond to 
physical differences in the velocities. The eigenvector 
(0, 1, —1) has eigenvalue 1/3 and corresponds to the per- 
turbation in which vp and v& are made unequal. So 
such a perturbation is irrelevant in the infrared, with 
scaling dimension — A* 2 /3 = —e/21. The final eigen- 
vector, (—12, 1, 1), induces a difference between the ve- 
locities of the fermions and bosons, and this perturba- 
tion is more strongly irrelevant: the scaling dimension is 
-13A* 2 /3 = -13e/21. 

So we have established that the Lorentz-invariant fixed 
point is at least linearly stable to velocity differences. 
However, the scaling dimension of the leading irrelevant 
operator is quite small (—e/21), and this could lead to 
very slowly decaying transients in the critical behavior. 



2. T > spectral functions 

The problem of computing T > 0, quantum critical, 
spectral functions in 2+1 dimensions is one of consider- 
able complexity. Eve n though the non-linear couplings at 
the fixed point, (|3.3|) , are small for small e, the low fre- 
quency response at T > cannot be computed in a bare 
e expansion. This is because the limits e — > and oj — > 
do not commute, and there are spurious low- frequency 
singularities in the bare e expansion: this is discussed 
further in Appendix [X| For the case of a scalar field 
alone, this problem of T > dynamics was addressed 
in Refs. [75],[l3]by a two-step procedure involving a map- 
ping to a quasi-classical problem. Here we shall use the 
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insight gained from these previous, controlled studies to 
motivate a simple self-consistent one-loop theory for com- 
puting the low frequency relaxation rates. In contrast to 
the previous studiesEj, we are aided here by the fact that 
damping appears already at the one- loop level in S s *+id 
or Sd 2 2+idx H i an d this simplifies the solution of the 
self-consistent equations. 

We shall be interested in the following-j-etarded re- 
sponse functions in the "quantum critical" £3 regime: 



G f (k,iw n ) = ZJ 1 ^i(fc,w n )*J(fe,w n )) 
G b (k,iuj n ) = Z^ 1 (4>(k,co n )(t>(-k, -0J n )) . 



(3.16) 



After analytic continuation to real frequencies Gf satis- 
fies the scaling form in (1.1), and an analogous result for 
Gb with / —y b. 



For w > T or fc > T, the response functions in (3.16) 



equal their T = results, and these can be computed in 
the bare e expansion. A familiar computation in the one 
loop renormalized theory gives us 



G f {k,u>) = Cfii- 



to + k x r z 



ky TX 



G b (k,Lu)=C b ^' n 



[k 2 - (w + zO) 2 ] 1 ~' ?//2 
1 



[k 2 - (u J + iO) 2 } 1 ' r,b/2 
where Cf t b are constants given by 

C f = 1 - A* 2 /4 + . . . 
C b = 1 - 2A* 2 + . . . . 



(3.17) 



(3.18) 



Clearly, (3.17) obeys the scaling form (1.1). 

Next we consider the more difficult case where u>, k 
are of order T or smaller. Here we shall proceed by us- 
ing an ansatz for the low frequency form of the Green's 
functions: this ansatz is motivated by the results of the 
corresponding low frequency regime in a number of other 
modelsLJ'Ej, especially the exactly solvable Ising chain in 
a transverse fielcO. In these studies, it was found that, 
provided the number of order parameter components was 
not large, the low frequency spectral functions had a sim- 
ple relaxational form, that of degrees of freedom in a 
simple dissipative medium. In the present situation, we 
have a single scalar order parameter, and can expect that 
a similar situation should apply. We therefore make the 
ansatz 

nti \ ( ^\^ s u + iTf + k x r z + k y T x 
G ^) = {t) B- {uj + l T f ) 2 



G 



i 



k 2 - (u + iT b )< 



(3.19) 



where T^f, are damping frequencies we need to deter- 
mine. The functional form for Gb in ( 3.19| ) was found 
to be a remarkably good fit to the low frequency portion 
of the known exact result for the, spectral function of an 
Ising chain in a transverse fieldOj'Ea. A further rationale 



behind (3.19) is that we have only included terms in the 
self energy of the fields which are lower powers of the fre- 
quency than those already present in the free propagator, 
and which are consistency with the requirement that all 
spectral functions are smooth at u> = for T > 0. It is 
precisely these low frequency powers which acquire a sin- 
gular form in the bare e expansion, and so we are forced 
to use the se lf-con sistent app roach to be described below. 
Comparing ( 3.19| ) with ( |l.l[ ) we see that consistency de- 
mands that be universal numbers times temperature. 
The purpose of the remainder of this subsection is to de- 
termine the universal numbers T^/T and Tb/T. Before 



turning to this, let us also quote, the form of (3.19) in 
imaginary frequencies 

/ /i \ -Vf i(j + iTfSga(oj n ) + k x r z + k y r x 
Gf(k,tuj n ) = (- 



J n\ +1/, 



G 



b (k,iuj n ) = (Jpj 



-m 



k 2 + (\Lu n \ + T b y 



(3.20) 



The damping represented by Tfb arises from interac- 
tions between the thermally excited bosons and fermions. 
The typical excitation will have energy of order T, and 
as the damping is dominated by the lowest energy exci- 
tations, the typical interaction vertex will have external 
frequencies at order T or lower. Motivated by this, we 
develop a perturbation theory for the interactions not in 
terms of the bare vertices, but in terms of the full inter- 
action vertices between the excitations. At the one-loop 
level, the damping terms arise from the Ao interactions 
alone: however rather than expanding in powers of Ao, 
we will express the self energy in powers of the full three- 
point irreducible vertex, A3 between one boson (</>) and 
two Fermi (^\ 2 , ^1,2) fields. A convenient choice is to 
use the zero momentum vertex between a Bose field at 
zero frequency, and Fermi fields at the minimum Mat- 
subara frequency of e n = irT. Bare perturbation theory 
for this vertex gives 



A3 — Aq — Aq 



d d k 
(2^p 



1 



(fc 2 + W 2)(fc 2 + (e„+ W „) 2 ) 
(3.21) 



where u> n is a fermionic frequency. [General scaling argu- 
ments from renormalization theoryEaLU show that at the 
fixed point the exact result for the vertex obeys 



A, 



T e/2 



z f z b 



T72 \ f) 



c. 



(3.22) 



where C3 is a universal number. Evaluatin g th e fre- 
quency sums and the momentum integrals in (3.21), and 
expressing everything in terms of renormalized couplings 
at the fixed point , and collecting low order terms in e, 
we find that ( 3.22 ) is indeed obeyed, with the universal 
number 
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Ca = 



S. 



1/2 
d+1 



(A* + 0.374367A 



*3 



(3.23) 



To proceed with the damping calculation, we write 
down the structure of the self-consistent one-loop equa- 
tions i n ter ms of the Green's functions of the bare fields; 
from (|le}) these are Gf = Z f G f and Gf = Z b G b . We 
define the self energies by 



A simple dimensional analysis of (3.27) shows that T can 
be completely scaled out of both equations for all d. The 
strength of the damping is determined by the dimen- 
sionless ratios Tf/T and T b /T, and these are completely 
determined b y the dimcnsionlcss universal C3. 

We solved ( 3.27 ) numerically in d — 2. The results are 
Tf/T = 0.581 and T b /T = 0.170. 



(Pf) 
(Gi) 



= -iu! n + k x r z + k y r x 
= ul + k 2 + s - Ef . 



yB 



(3.24) 



The self-consistent, one-loop expression for the self ener- 
gies, expressed in terms of A 3 , are 



E?(fc,w n ) 



X Gf (k-PiUJr, 
1(1 , 



en) 



e n ) (3.25) 



(2tt) q 
x Gf(k +p,uj n 



We will now express this in terms of renormalized quanti- 
ties. First we note by comparing (3.24) with (3.19,3.20) 
that we only need the imaginary part of the self ener- 
gies at small real frequencies; in particular, the damping 
co-efficients can be expressed as 



Zt lim ImS? (0,w) 



r b =(»Y Vb Z b lnn Im ^ ^ 
\TJ w-Q 2w 



(3.26) 



We now insert ( ^20| ), ( gTgg ), d3.25[ ) into ( |3~26| ). To our 

satisfaction, we find that all factors of the renormaliza- 
tion factors Z b j and the scale [i precisely cancel out, and 
the remaining expressions involve only univers al qu anti- 
ties. Perfor ming the frequency summation in (3.25) and 
inserting in ( p. 26] ) we obtain 



3 ' 7T J (2ir) d sinh(n/T) 
1 



ci — 

J I 

xlm 
xlm 



dtt f d d k 



T< 



2C\ 



k 2 - (fi + iT b ) 2 

n + iTf 

k 2 - (fl + iTf) 2 

00 dn r d d k r e_1 

7T 



Im 



(27r) d cosh 2 (f7/2T) 

2 



n + ir 



f 



-k 2 I Im 



k 2 - {n + iT f ) 2 

2 ' 



1 



k 2 - (Sl + iTf)' 



(3.27) 



B. C symmetry breaking in a d-wave superconductor 

A number of transitions involving C symmetry break- 
ing in a d-wave superconductor were noted in Section ||. 
These involve the onset of either stripe or spin-Peierls or- 
der, and such transitions appear in Figures ||, fj], and [n]. 
In all cases, the order parameter can be identified with 
scalars $ x , <& y representing the amplitude of charge den- 
sity waves with wavevectors (Q,0) and (0,(2). If <2 is 
commensurate with the underlying lattice, then & x ,y are 
real; otherwise & x ,y are complex, with their phases rep- 
resenting the freedom of the charge density wave to slide 
with respect to the underlying lattice. On general sym- 
metry grounds, we c an write down an effective action for 



<& x ,y, similar to (3.3): 



S 1 * = / d d xdr 



\d T 3> x \ 2 + \d T 3> y \ 2 + |V$a| 2 + |V$i 



-so(|<M 2 + |^| 2 ) + ^(|^| 4 + |$,| 4 ) 



(3.28) 



ir X)2/ can be either real or complex. First- 



where, for now, <E> 2 
order time derivative terms, like Q x d T $ x , are forbidden 
here by spatial inversion symmetry under which <fr x — > 
and such a term changes sign after integration by 
parts. 

To complete the theory, we have to consider the cou- 
pling of *&x,y to the gapless Fermi excitations at wavevec- 
tors (± K, ± K) . Conservation of momentum implies that 
there is in fact no long-wavelength coupling between & x . y 
and (which is linear in $ x ,y) unless Q = 2K . The 
mean-field studies of Section [n] always obtained Q ^ 2K, 
and this expected to be the generic behavior. How- 
ever, we cannot rule out the possibility that there is 
a mode-locking phenomenon which preferentially con- 
denses a charge density wave at wavevector Q = 2K over 
a finite range of parameters. 

For Q 7^ 2if, 5$ is the complete critical theory of the 
transition: the fermions are not part of the critical theory 
and so the transition is in class B. The simplest allowed 
couplings between the fermions and the critical degrees 
of freedom are terms like wq J d d xdT\Q x \ 2l i>\T z '$>i. Sim- 
ple power counting shows that wq has scaling dimension 
1/u — d, where v is the c orrela tion length exponent of the 
transition described by (3.28). We expect that this v is 
greater than that of the rf+1 = 3 dimensional XY model, 
which is rs 2/3, and hence Wq is irrelevant in d = 2. By 
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counting scaling dimensions (or by an explicit perturba- 
tive computation) we can deduce that the self energy of 
the nodal fermions obeys ImE^ ~ w\T 2d+1 ~ 2 l v , and so 
the damping rate vanishes with a super-linear power of 
T as T — ► 0, as expected for a class B transition. 

In the remainder of this section, we consider the class 
A transition with Q = 2K, and both incommensurate 
(experimentally, and in our mean-field theory, K is in- 
commensurate), so that $> x .y are complex. Now a cou- 
pling between <& x ,y and ^1.2 is possible. Writing dowa 
all possible terms consistent with symmetries we obtainEl 

SW = Jd d xdT (Ao+Co) {^>x^2 TZ ^l+^y^ab^2aT X ^ 
- (Ao ~ Co) (Ktfjr"*! + <S>yE ab ^2aT Z ^U 



H.c. 



(3.29) 



The renormalization group analysi s of + 5$ + 
parallels that carried out in Section III A . and so we will 



be brief. The theory can be shown to be Lorentz invari- 
ant for c — vp — va and Co = 0, and so also has z = 1. 
For this Lorentz-invariant case, the renormalization con- 
stants, replacing those in ( 3 . 1 0| ) are 



z b = 


1 - 


2A 2 
e 




Zf = 


1 - 


A 2 
e 




Zx = 


1 - 


A 2 
e 




Z u = 


1 + 


5u 2 - 


h v 2 - 2A 4 




uc 


Z v = 


1 + 


2v 2 h 


- Auv - 4A 4 



ve 



(3.30) 



(the coupling v has been define d fro m v$ following the re- 
lationship between u and uq in (3.9)), the beta functions, 
replacing those in (3.11) are 



/3(A) 



-A + A 



(3(u) = -eu + 5u 2 + v 2 + 4uA 2 - 2A 4 

0(v) = -ev + 2v 2 + Auv + Av\ 2 - 4A 4 , (3.31) 



and the anomalous dimensions modifying those in ( 3.13| ) 



are 



rib = 2\* 2 

\ *2 

ri f =\ . 



(3.32) 



(Hi) In (3.27), the integrand in the expression for Tf 
has a prefactor 2C|, while that for has a prefactor 
Cf . The numerical values of the damping coefficients are 
now Tf/T = 1.35 and T b /T = 0.395. 



C. OAF order in a d-wave superconductor 

We discuss here the transition between a d- wave super- 
conductor and the phase with coexisting OAE.Iar stag- 
gered flux) and d-wave superconducting ordert3'E3 : such 
a transition appears in Fig. |ll| at S w 0.12. We will 
show th at t his transition is in class B (in the notation of 
Section LA). r-. 

As has been emphasized by NayakcJ, the OAF is char- 
acterized by a d-w&ve order parameter in the particle-hole 
channel: 



J 



Ck.a) 



(cos k x — cos k y 



(3.33) 



where G = (jr, n) and cj> is a real order parameter. As 
in the above subsections, the key issue is the coupling 
of this order parameter to the fermionic quasiparticles 
of the (i-wave superconductor. As the order parameter 
carries momentum G, no coupling, linear in <f> is possible 
unless the nodal points are at (±K,±K) with K = ir/2. 
For K 7^ 7r/2 the simplest allowed coupling is 

w J d d xdT(j} 2 ^\T" 



+ As we saw in Section IIIB , 

is irrelevant and places the transi- 
The damping rate obeys ImS^ ~ 



such a coupling 
tion in class B. 
w 2 T 2d+i-2/^ _ where Vj _ 63 ig the corre i a . 

tion length exponent of the Ising mode l in d + 1 = 3 
(described by the field theory in (3.3)). 

For completeness, let us also cons ider the special case, 
where K = n/2. Then, from (3.33) we can computet^ 
the following coupling between <fi and ^1,2: 



0£ab 



ity la T x d v ty lb + i^2aT X 3 x ^2b + H.c. 



(3.34) 



Note that this coupling has one more derivative than 



those in (3.6) and (3.29); this is a conseq uence of the 
vanishing of the factor (cos k x — cosk y ) in (3.33) at the 
nodal points. Therefore, (3.34) is ir relev ant by simple 



power- counting. A coupling such as (3.34) will not lead 
to a fermion spectral function obeying (1.1); instead the 
imaginary part of the fermionic self-energy vanishes as 
ImE^ ~ r p 2 +'m 5 where r\i > is anomalous dimension 
of the order parameter ((f)) of the Ising model in D = 3 
spacetime dimensions. So even for K = tt/2, this transi- 
tion remains in class B. 



The computation of the T > spectral functions pro- 
ceeds as be fore, but with the following changes: (i) In 
( |3.21 ) and (3~23| ), the A 3 terms have the opposite sign 



(ii) In ( 3.25 ), the integrand in the expression for has 



a prefactor 2A|, while that for E^ has a prefactor —A 



IV. DISCUSSION 

This paper has presented a comprehensive mean-field 
study of realistic models of the cuprate superconductors. 
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A representative sample of our results appears in Figs ||- 
[ll], and the properties of the phases therein were summa- 
rized in Section |. Xh csp ^aj^dd-yesults. unify many 
other earlier studiefiettaSBffiaffie and expose 
the relationships between them. 

A second focus of the paper has been on the second- 
order quantum phase transitions in Figs We 
paid particular attention to the T > fermionic quasi- 
particle spectra in the vicinity of the nodal points in 
the Brillouin zone, with the purpose of understanding 
the observed quaatum-criticality in recent photoemis- 
sion experiments^. We divided the transitions into two 
classes, A and B (described in Section I A), with only 
those in class A leading to universal damping with lo/T 
scaling near the nodal points. Particularly appealing ex- 
amples of class A transitions, for which class A behav- 
ior was generic and did not require any special param- 
eter values, were those involving time-reversal (T) sym- 
metry breaking in a d-wave superconductor: the most 
important of these are the transitions from d-wave to 
(s* + id)-wave or (d x 2_ y 2 + id xy )-wave superconductiv- 
ity, which appeared in our mean- field phase diagrams. 
The transition to (d x 2_ y 2 + id xy ) order had the addi- 
tional satisfying feature of very naturally leading to the 
absence of quantum-critical damping of quasiparticles at 
momenta (tt, fe,7r) (with < fc < n), as is found 
in experimentaljcZI below the superconducting critical 
temperature. However, the transition involving onset 
of "staggered- flux" (or orbital antiferromagnet) order in 
a d-wave superconductor, which broke both T and C 
symmetries, was not of class A. We note that fermion 
damping in a model involving T symmetsy breaking has 
also been examined recently by VarmaE3, although he 
refers to quantum criticality associated with a transition 
in a Fermi liquid and not a superconductor. We also 
examined the onset of "charge stripe-order" (C symme- 
try breaking) in the d-wave superconductor: such transi- 
tions belonged to class A if the charge ordering wavevec- 
tor was precisely equal to the separation between two 
nodal points of the d-wave superconductor. This is a 
fine-tuning condition, which is also not supported by ex- 
periments, and makes the C breaking transition a less 
attractive scenario for explaining fermion damping. 

So the most viable candidate for the state X in Fig || is 
the (d x 2_ y 2 + id^-wave superconductor. For this case, 
the damping mechanisms appear to divide the fermion 
excitations into two distinct components. The fermions 
along the (1,0), (0,1) axes are strongly paired in the 
d a .2_ !/ 2-wave state but are decoupled from the critical or- 
der parameter fluctuations (d xy ) to the state X: con- 
sequently there is negligible damping of these fermions 
below T c . On the other hand, the fermions along the 
(1,±1) axes couple strongly to the d xy order parame- 
ter and undergo quantum critical damping as described 
by ( |l . 1| ) . The situation changes dramatically once we 
go above T c . Now phase fluctuations and the prolifera- 
tion of hc/2e vorticesc2l will strongly scatter the fermions 
which couple efficiently to the predominant d x 2_ y 2 or- 



der parameter: these are the fermions along (1, 0), (0, 1) 
directions, while the vortices are largely invisible to the 
fermions along the (1,±1) directions. Moreover, the an- 
tiferromagnetic spin fluctuations, which were responsi- 
ble for fermion pairing along the (1, 0), (0, 1) axes below 
T c , will scatter these same fermions (on the "hot spots") 
above T c ; again, these fluctuations are invisible to the 
(1, ±1) direction fermions because the antiferromagnetic 
wavevector does not co nnec t the nodal points. Indeed, as 
we indicated in Section I A , the predominant damping of 
the (1, ±1) direction fermions above T c continues to arise 
from the quantum critical d xy fluctuations to the state X: 
this mechanism applies as long as the quantum-critical 
scattering length of these fermions remains shorter than 
the superconducting phase coherence length. 

With an eye, towards comparisons with photoemission 
experimentalj'Ha, we review our results for the nodal 
fermionic spectral functions of the class A transitions. 
As in Ref. [lC] we will follow the evolution of the spectral 
function along a line from the zone center going through 
the nodal points e.g. from (0,0) along the (1,1) direction 
through the nodal point at (K,K). At the wavevector, 
(K + k,K + k), our results for ( L.l) are contained in the 
diagonal components of ( [3.17 , 3 . 1 9| ) ; we express these 
results in the form 



-Vf 



(4.1) 



Note that E/ is strictly not a self-energy (and thus the 
tilde), as some of the self-energy corrections have already 
been abs orb ed into the prefactor of an anomalous power 
of T in (4J). We have also set the velocity in the (1,1) 
direction to unity. 



For small u>, k, our result for £/ in (3.19) is 



U and \k\ < T, 



(4.2) 



where is a universal number. For the d-wave to 

(s* + id)-wave or (d x 2_ y 2 + id X y)-wave transition we esti- 



mated V f/T — 0.58 in Section III A 2, while for the onset 
of a certain type of C ordering in a d-wav e sup erconduc- 
tor we obtained Jj/T = 1.35 in Section [HIB] . We plot 
the results (OO) in Fig. fil. 



For large u> or fc, our result is in ( 3. 17 ) . For small rjf, 
this can be written as 



2 



(k-u) [ In 



or 



j,2 



- 1 



+ i-k\lo - k\6{\u>\ - \k\ 



M or |fc| > T, 



(4.3) 



where 9 is the unit step function. Note that the imagi- 
nary part of vanishes for |u;| < \k\, and, in the present 
form, this will lead to an infinite spectral density at the 
threshold \lo\ = \k\. H owe ver, this is repaired by con- 
sidering corrections to ( |4.3| ). Within the scaling limit of 
the universal theor ies b eing considered here, we evalu- 
ate the expression (3.25) in Appendix [AJ at T > 0, but 
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FIG. 12. Low-frequency photo-emission intensity nea r 
a nodal point, Im Gf(k, uf) n/(o»), given by (p~l|), (4.2). 
Here, Tf/T = 0.58, which is the result for the d-wave to 
(s* + id)-w ave or 
Sec. 



_ y 2 + id xy )-wa.ve transition discussed in 
III A, and nf(u) = [exp(tj/T) + l]" 1 denotes the Fermi 



function. 



with 



Ifel ^> T, we obtain in addition to the lead- 



ing term in (4.3), subleading T-dependent corrections 
ReE^ - T 2 /uj and ImE^ - T 3 /lj 2 . These are still 
very small contributions, and so we can expect that the 
system will be exceptionally sensitive to non-universal 
corrections to scaling right at the threshold frequency. 
We believe that the most important of these will come 
from elastic scattering off impurities; for a weak impu- 
rity scattering potential, U lmpi we have the additional 
contribution 



^impM 



(4.4) 



We have added (4.4), with a very small prefactor, to (4.5) 
and plotted the result in Fig. [l3|; the contribution of (4.4) 
can be neglected almost everywhere except right near the 
threshold. 

It is interesting to note that our results (4.2) and (|4 



J) 

bear a superficial similarity to the "marginal Fermi liq- 
uid" fitting functionslEJ. More specifically, (i) in the lat- 
ter approach, the prefactor of the power of T in fl4.l| ) is 
absent; (ii) the small |w|/T behavior of the self-energy 
in (4.2) is similar to that of the marginal Fermi lin - 
uid; (Hi) for large \u)\/T, the k dependence in ( |4.3| ) 
(as in ImE ~ \ui — k\6(\u>\ — \k\)) is replaced simply by 
ImE ~ \uj\. A significance consequence of this last dif- 
ference is that our large k spectral densities have more 
asymmetric lineshapes (even before being multiplied by 
the Fermi function, as is necessary for photoemission ex- 
periments), than those found in the marginal Fermi liq- 
uid functions; this is illustrated in Fig. 



14 



It would be 

interesting fo r ex periments to test for the k dependence 
predicted in (4.3). 




a/T 



FIG. 13. High-frequency photoemissio n in tensity near a 
nodal point, 1xaGf{k,uf)nf(uf), given by ( |4.3| ), (4.4), with a 
small impurity contribution L^ np =0.1. 
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FIG. 14. Sample spectral function for k/T = 10 - this 
is a cut through the spectrum of Fig. |l^. The asymmetry of 
the lineshape is clearly visible (and it is not simply due to the 
Fermi function occurring as prefactor of the photoemission 
intensity) . 
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Finally, we reiterate an important feature of our anal- 
ysis of quantum criticality: we have only found Lorentz- 
invariant fixed points at which a ll excit ations have equal 
and isotropic velocities. Section III A 1 showed that, for 
the d-wave to (s* + id)-w&ve or (d x i_ y i + id xy )-w&ve 
transition, such a fixed point was at least linearly sta- 
ble to perturbations which e.g. set vf ^ va- However, 
this does not rule out the possibility that there may be 
other non-Lorentz-invariant fixed points of the renormal- 
ization group equations in which vpb^A is significantly 
different from unity. Experimentallyea, it is clear that 
vf va, but this is still compatible with a Lorentz- 
invariant fixed point: we found in Section [II A 1 that the 



leading irrelevant operator which breaks Lorcntz symme- 
try had a scaling dimension of very small absolute value 
(e/21 w 0.048), and so the system can reside in a tran- 
sient region with vf over a very wide temperature 
range. This issue will be addressed further in future work. 

Notes added: (i) Recent THz conductivity measurements 
on Bi 2 Sr 2 CaCu2 8 +5 by Corson et allEi have obtained a 
quasiparticle relaxation rate linearly proportional to T, 
at temperatures well below T c . Combined with the pho- 
toemission experimentsCJ, these results provide strong 
support for quantum critical damping of the nodal quasi- 
particles in the d-wave superconductor, as is expected 
near a class A quantum critical point between two su- 
perconducting states. Corson et al. observe the quantum 
critical damping above 20K, which suggests that the en- 
ergy per coherence volume of the second superconducting 
state (say, the (d x 2_ y 2 + id xy ) superconductor) is higher 
than that of the <i-wave superconductor by less than 20 

(ii) A recent workU has given a unified discussion of the 
quantum phase transitions considered here, and those in- 
volving M. symmetry breaking considered in Ref. p9[ 
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APPENDIX A: PERTURBATION THEORY FOR 
FERMION DAMPING 

This appendix will compute the T > scaling function 
for the fermion spectral weight in (1.1) for the case of a 
quantum critical point between a ci-wave and a (s* + id)- 
wave or (d x ?_ y i + id xy )-w&ve superconductor. We will 
use a simple renormalized perturbative expansion at the 
fixed point found in Section III A. We will show that 



such a procedure leads to spurious singularities in the 
low frequency regime htu < ksT. These singularities 
were cu red by the self-consistent analysis described in 
Section til A 2 . 



From the expression (3.25), we obtain to leading order 
in A 2 and e 



-Vf 



G / 1 (fc,w„) = 
{-iuj n + k x r z + k y T x ) 



1*2 



+1 



(2n) d ^ 



—it. 



PxT 



ln(/VT)J 



(p 2 + 4)((fc-p) 2 + K-e«) 2 )' 



Evaluating the frequency summation, and performing the 
momentum integration over the terms not involving any 
thermal Bose or Fermi factors, we obtain 



(|) V! Gj 1 (k,u n ) = -zf\k,w n ) + 



{-iw n + k x r z + k y r x ) I 1 - 2£ I In 



eT 2 



(Al) 



where we have used ( 3.13 ), and is a thermal contri- 
bution which vanishes as T — * 0. To leading order in e, 
the expression for si 1 can be evaluated in d = 3, and 
we obtain for k along the x direction (recall that below 
( |3 . ip we rotated the axes by 45 degress from the axes of 
the square lattice): 



tf{k,u n ) = -16^77/ 



d 3 p 



1 



WlP- 



(2tt) 3 2p|p - k| 
k|) + /(p)][(p»r*(p-|p-k|)-»w n p] 



+ 



(p-|p-k|) 2 
[n(\p--k\)-f(p)}[ Px T z (p+\p- 



(P- 



■U3* 



(A2) 



where p = |p|, k = |k|, n(k) is the Bose functio n, a nd 
f(k) is the Fermi function. The expression in (A2) is 
reliable for \u>\,k ^S> T: evaluating the integra ls in this 
regime we obtain the estimates quoted below (4.3). On 
the contrary, for \u>\,k <^ T, the above expressions are 

pathological; we obtain e.g. ImT,f\o,Lj) - T 2 5{u). This 
should be contrasted with t he smooth behavior as a func- 
tion of lu assumed in ( |3. | LiJ ). The latter is the correct re- 
sult on physical groundsH, and estimation of the damp- 
ing constants requires a self-consistent approach like that 
followed in Section III A 2| . 
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